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ABSTRACT 


This  report  presents  a  summation  of  the  work  performed  under  the 
Antenna  HF  Optimization  Study,  a  study  to  determine  the  optimum  log  periodic 
dipole  ar^ay  design  for  a  high  gain,  low  coverage  HF  antenna  over  real  and/or 
ground  screen  covered  earth.  The  major  problem  encountered  is  one  of  deter¬ 
mining  antenna  electrical  characteristics  accurately  in  free  space  and  in  a  real 
earth  environment.  The  present  state  of  the  art  in  analytical  antenna  design 
does  not  allow  the  engineer  to  predict  log  periodic  antenna  electrical  character¬ 
istics,  especially  over  real  earth,  with  sufficient  accuracy  to  be  valuable  for 
strict  system  requirements.  To  solve  this  problem,  a  mathematical  model  of 
a  log  periodic  ante)_ia  is  set  up  and  solved  with  the  use  of  high  speed  digital 
computers.  The  model  used  allows  one  to  obtain  the  current  distribution  on 
the  antenna  knowing  only  its  physical  dimensions  and  something  about  the  con¬ 
ductivity  of  the  material  used  in  its  construction.  The  theoretical  model,  which 
treats  the  antenna  as  a  boundary  value  problem,  automatically  includes  all  the 
interactions  or  mutual  coupling  between  elements.  From  the  computed  values 
of  the  currents  on  the  antenna  all  electrical  characteristics  of  the  antenna,  such 
as  radiation  patterns,  absolute  gain,  and  input  impedance  are  calculated. 
Theoretical  and  measured  UHF  radiation  patterns  and  gains  are  compared  for 
evaluation  of  the  technique.  Free  space  and  real  earth  calculations  are  presented 
and  an  optimum  antenna  design  is  chosen  from  this  information. 
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EVALUATION 


1.  The  objective  of  this  effort  was  no  investigate  antenna  techniques 
applicable  to  arrays  of  HF  dipole  let  periodic  elements  to  provide  gains 
of  17  dbi  and  radiation  coverage  with  a  minimum  of  14  dbi  from  two  (2) 
through  16  degrees  in  elevation  and  plus  and  minus  8  degrees  in  azimuth 
over  the  3-30  MH2  frequency  range  with  minimum  real  estate. 

2.  The  approach  taken  by  Sylvania  was  to  investigate  H-Plane  arrays  of 
log  periodic  elements  of  the  DuHamel  and  Mei  configurations.  The 
original  approach  was  to  perform  both  an  analytical  and  experimental 
effort  to  determine  the  advantages  of  the  Mei  array  approach  over  the 
conventional  DuHamel  array  configuration.  After  the  effort  had  been 
initiated,  it  was  redirected  to  put  full  emphasis  on  developing  a 
computer  program  which  would  have  the  capability  of  computing  the 
performance  characteristics  of  dipole  log  periodic  elements  and  arrays 
of  log  periodic  elements  in  both  the  DuHamel  and  Mei  configurations. 

This  computer  program  was  written  for  the  RADC  GE  635/64S  Computer  in 
order  that  RADC  antenna  engineers  could  make  efficient  use  of  the  computer 
program  using  RADC's  computer  facilities.  RADC  used  its  computer 
facilities  with  the  Sylvania  antenna  program  to  obtain  the  performance 
characteristics  of  log  periodic  arrays  as  presented  in  this  report.  It 
was  determined  by  the  many  computer  runs  performed  by  RADC,  that  the 
computer  program  developed  by  Sylvania  has  a  limitation  of  only  being 
able  to  handle  log  periodic  structures  with  bandwidths  less  than  an 
octave.  The  results  from  analyzing  lor  periodic  structures  larger  than 

an  octave  bandwidth  did  not  correlate  very  well  with  experimental  results. 

3.  The  results  of  this  effort  have  shown  that  the  following  parameters 
have  a  significant  effect  on  the  performance  of  dipole  log  periodic 
elements: 


a.  Dipole  diameter 

b.  Characteristic  impedance  of  feed  line 

c.  Tau  factor 

d.  Alpha  factor 

e.  Number  of  dipo’e  elements  on  structure 

f.  Bandwidth  of  the  structure 


ix 


This  co^mter  program  will  enable  the  antenna  engineer  to  determine 
optimum  antenna  parameters  for  arrays  of  dipole  log  periodic  elements 
before  an  extensive  experimental  effort  is  undertaken.  It  is  hoped 
that  the  limitations  of  this  program  can  be  overcome  and  can  be  made 
to  satisfactorily  handle  greater  bandwidth  structures. 

4.  With  the  use  of  the  RADC  GE  635/645  Computer  Facilities  and  the 
computer  program  developed  under  this  effort,  the  Air  Force  will  be 
able  to  analytically  determine  the  performance  characteristics  of 
large  HF  log  periodic  dipole  arrays  for  future  radar  systems.  This 
capability  will  also  be  advantageous  for  analyzing  in  more  detail 
the  performance  characteristics  of  dipole  log  periodic  structures. 


CARMEN  S.  MALAGISI 
Project  Engineer 
RADC/EMATA 
Griff iss  AFB,  NY 


I  INTRODUCTION 


Insufficient  information  on  the  electrical  characteristics  of  antennas 
is  and  has  been  a  continuous  problem  to  the  systems  engineer.  Whenever  he 
attempts  to  predict  the  receiver  signal  level  a  lack  of  detailed  gain  character¬ 
istics  becomes  significant.  The  amount  of  information  available  varies  greatly 
from  one  type  of  antenna  to  another.  The  antenna  we  wish  to  consider  in  this 
study  is  the  log-periodic  antenna. 

There  are  may  different  forms  of  log-periodic  antennas;  however, 
only  one  may  be  considered  as  useful  at  the  present  time  for  the  operating 
frequencies  many  fo  r  long-range  purposes  --  the  log-periodic  dipole  antenna 
(L.PDA).  Figure  1-1  is  an  example  of  a  VHF  log-periodic  dipole  antenna,  the 
absolute  gain  of  which  is  too  low  for  long-range  systems.  These  antennas  may 
of  course  be  arrayed  to  increase  the  gain  to  the  required  level.  Arraying  may 
be  accomplished  such  that  the  resulting  array  has  the  same  frequency-indepen¬ 
dent  characteristics  as  the  fundamental  antenna.  This  technique  has  a  theoretical 
maximum  gain.  Figure  1-2  is  an  artists  concept  of  such  an  array  in  the  HF 
region.  One  may  also  array  the  antennas  to  give  broad-band  characteristics 
that  change  with  frequency,  but  have  no  theoretical  limit  on  maximum  gain; 
however,  at  the  present  time  we  are  concerned  only  with  the  frequency-indepen¬ 
dent  array. 

The  determination  of  the  electrical  characteristics  of  a  LPDA,  when 
the  antenna  is  a  close-spaced  array  on  a  homogeneous  earth  covered  by  a 
ground  screen  of  finite  length,  is  not  an  easy  task.  The  first  step  in  the  complete 
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Figure  1-1.  VHF  Log-Periodic  Dipole  Antenna. 


L _ (Continued) 

solution  of  this  problem  is  the  determination  of  the  theoretical  characteristics 
of  such  an  antenna  in  free  space.  This  report  is  addressed  primarily  to  this 
problem.  The  mathematical  techniques  to  be  developed  must,  of  course,  be 
of  a  form  that  is  applicable  when  the  complete  antenna  is  composed  of  arrays 
of  LPDA's. 

The  LPDA  is  in  a  class  of  antennas  that  to  the  present  day  do  not  have 
suitable  approximate  mathematical  solutions.  Attempts  to  compensate  for  the 
mutual  coupling  between  the  non -equal -length  close-spaced  dipole  elements 
using  the  usual  arraying  theory  have  been  unsatisfactory.  The  present  state 
of  the  art  in  analytical  antenna  design  does  not  allow  the  engineer  to  truly 
predict  the  performance  of  log-periodic  dipole  antennas. 

The  assumptions  of  the  form  of  the  current  distribution  on  the  dipole 
elements  and  consideration  of  the  interaction  between  pairs  of  dipoles,  neglecting 
the  effects  of  neighboring  dipoles,  has  given  useful  but  limited  results. 

The  techniques  used  in  obtaining  solutions  of  the  LPDA  go  beyond  the 
usual  analytical  methods  and  make  maximum  use  of  advanced  numerical 
analysis  and  high-speed  digital  computers.  Many  "solutions"  to  antenna  prob¬ 
lems  in  the  past  have  made  assumptions  on  the  current  distribution  on  the  elements 
and  have  used  this  assumed  current  distribution  to  determine  radiation  patterns, 
etc.  In  other  -vords,  an  answer  has  been  assumed  for  a  major  portion  of  the 
problem  and  only  the  minor  problem  of  determining  the  far -field  pattern  (given 
the  current  distribution)  has  been  solved.  In  this  report  no  assumptions  on  the 
distribution  of  current  are  made;  the  antenna  is  treated  as  a  boundary -value 
problem  for  which  a  solution  can  be  obtained  by  numerical  analysis  and  the  aid 


of  a  digital  computer. 


L _ (Continued) 

Thus  only  those  parameters  that  define  the  antenna  physically  and 
the  operating  frequency  are  the  input  parameters  to  the  digital  computer  pro¬ 
gram.  The  computer  program,  in  turn,  determines  the  current  distribution 
on  the  antenna  based  on  these  physical  parameters.  From  these  derived  values 
of  the  currents  on  the  antenna,  all  electrical  characteristics  of  the  antenna, 
such  as  radiation  patterns,  maximum  absolute  gain,  and  input  impedance  are 
calculated. 

As  shown  in  Section  III  the  boundary-value  problem,  which  mathematically 
represents  the  antenna,  can  be  represented  by  a  Fredholm  integral  equation. 

(A  Fredholm  integral  equation  is  one  where  the  unknown  is  under  the  integral 
sign. )  It  can  be  shown  that  such  an  integral  equation  can  be  considered  as  the 
limit  of  a  set  of  n  algebraic  equations  where  the  accuracy  of  the  solution  is  a 
function  of  n.  Use  is  made  of  this  fact  by  writing  the  n  algebraic  equations 
and  solving  the  n  equations  using  matrix  techniques.  These  equations  are, 
of  course,  complex  in  the  mathematical  sense. 

Th'-  next  section  of  this  report  presents  the  theory  behind  the 
mathematical  model  of  an  n  element  log-periodic  dipole  antenna  in  free 
space.  The  third  section  presents  some  of  the  numerical  results  obtained 
from  the  computer  program  which  utilized  the  mentioned  mathematical  model. 
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II  THEORY 


A. _ INTRODUCTION 

The  primary  objective  of  this  section  is  to  present  the  theory 
behind  th:  mathematical  model  of  an  n  element  log-periodic  antenna  in 
free  space.  The  analysis  of  this  mathematical  model  by  numerical  tech¬ 
niques  will  be  discussed  simultaneously  to  enable  the  reader  to  understand 
the  capabilities  and  limitations  of  the  computer  programs  utilizing  this 
particular  model.  This  section  does  not  provide  a  documentation  of  the  pro¬ 
gram  that  has  been  written  using  the  model  to  be  discussed,  but  merely 
provides  the  theory  and  methods  used  to  write  such  a  program. 

The  next  part  of  this  section  provides  a  brief  outline  of  the  basic 
theory  involved  in  the  mathematical  modeling  of  a  wire  antenna  beginning  with 
a  simple  dipole  and  then  extending  this  concept  to  include  a  log-periodic 
antenna.  The  remainder  of  this  section  provides  a  detailed  analysis  of  how 
the  model  is  derived  from  the  theory  and  what  numerical  techniques  are  used 
to  implement  this  theory  on  a  digital  computer. 
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B.  OUTLINE 

To  find  an  expression  for  the  current  on  a  dipole  in  terms  of  the  applied  voltage  it  is 
convenient  to  first  derive  an  equation  relating  the  voltage  to  the  vector  potential,  and  then 
find  an  expression  for  the  vector  potential  in  terms  of  the  current. 

The  electric  field  at  any  point  outside  a  dipole  can  be  expressed  as 


E  =  -  V0  -  jooA  . 


Now,  defining 


v.A  axv 


V-A  -  -jcopcd  , 


(1) 


(2) 


we  get 


E  = 


i-  [„ 

jpe  l  ' 


jcope 


V-A  +  «2  pe A)  ] 


(3) 


Assuming  the  dipole  parallel  to  the  z-axis.  Equation  (3)  reduces  to 


d2A  „ 


(4) 


where 


.2  2 

k  =  ui  Pc  • 


For  a  perfect  conductor  the  tangential  electric  field  is  zero  on  the  boundary.  Equation  (4) 
then  becomes 


A\  2 
— #  +  k2Az  =  0 


(•r») 


dz 
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B  ■  (Continued) 


The  general  solution  of  this  differential  equation  is 

A  =  C.  cos  kz  +  C-  sin  kz  .  (0) 

z  1  2 

With  a  slice  generator  of  voltage  V  applied  at  the  center  of  the  dipole,  as  derived  by 

/  5 

Hallen,  Equation  '<!)  becomes 


V  I  ! 

Az  =  C i  cos  kz  -g-  sin  kjz|  (7) 

This  is  the  vector  potential-voltage  relationship  that  we  set  out  to  find. 

To  relate  the  voltage  to  the  current  flowing  on  the  dipole,  we  can  now  use  the  known 
relationship  between  the  vector  potential,  A  ,  and  the  dipole  current, 


A 


z 


J£ 

4tt 


e-jkR 

R 


dz’  , 


(8) 


where  R  is  the  distance  between  the  point  in  space  and  the  point  on  the  dipole  where 

current  I«  is  flowing, 
z 

Combining  Equations  (7)  and  (8)  we  obtain  the  desired  current-voltage  relationship. 


g-jkR  ,  y 

— — —  dz  =  B  cos  kz  +  j-^sinkjzl 


(9) 


In  order  to  solve  this  integral,  we  approximate  the  integral  by  a  finite  sum  of  currents 
along  the  dipole.  Each  term  in  this  sum  contains  an  unknown  value  of  current.  To  solve 
for  these  unknown  current  values,  we  need  the  same  number  of  equations  as  we  have  incre¬ 
ments  of  currents  along  the  dipole.  By  specifying  the  value  of  z  to  a  specific  increment 
along  the  dipole,  one  equation  can  be  written  for  each  current  segment.  With  the  dipole 
divided  into  N  segments,  we  now  have  N  number  of  independent  equations .  Atypical 
equation  for  the  i-th  segment  is 


H-3 


B.  (Continued) 


where  n  covers  all  dipole  segments,  including  the  i-th  segment.  In  addition  to  the  N 
number  of  unknown  currents,  the  integration  constant  B  is  unknown. 

An  additional  independent  equation  is  obtained  from  the  boundary  condition  that  the 
current  at  the  tip  of  the  dipole  is  zero; 

■  °  •  (“i 

With  N4 1  independent  equations,  the  N-+1  unknowns  can  now  be  determined. 

The  same  theory  can  be  applied  to  a  series  of  dipoles ,  as  in  the  case  of  a  log-periodic 
dipole  array.  For  a  log-periodic  array  with  M  number  of  dipoles.  Equation  (10)  becomes 


M  N 

E  E  i2 

m=l  n=T 


t  V-  I 

—  bz  =  B.  cos  kz..  +  j-^  sink  z..  , 

mn  j  )i  3  2.7)  I  Ji 


where  m  covers  all  dipoles,  including  the  j-th  dipole. 

The  unknown  voltage,  V.  ,  at  the  base  of  the  j-th  dipole  in  the  array  is  related  to  the 
applied  voltage,  V  ,  at  the  input  terminal  of  the  array  by  a  series  of  tranmission  line 
equations . 

The  currents  on  the  dipoles  of  a  log-periodic  antenna,  calculated  by  using  the  approach 
described  above,  are  by  definition  the  actual  currents,  including  all  mutual  interactions 
between  the  dipoles. 

From  these  values  of  currents  on  the  structure  all  electrical  characteristics  of  the 
antenna,  such  as  far-field  patterns,  the  absolute  gain  and  input  impedance  can  be 
calculated. 
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C. _ THE  DIPOLE  ANTENNA 

The  cylindrical,  cer  '-.er-driven  antenna  has  been  analyzed  as  a 
boundary  value  problem  by  many  workers.  Basically,  two  methods  of  attacking 
the  problem  of  determining  the  distribution  of  current  along  the  center-driven 
antenna  have  been  used.  In  the  first  of  these,  which  is  not  a  boundary  value 
treatment,  one  uses  the  apparent  similarity  between  the  antenna  and  an  open- 
end  parallel- wire  transmission  line.  Thus,  by  suitably  correcting  existing 
transmission  line  theory,  a  satisfactory  approximation  for  the  antenna  can  be 
determined. 

The  second  method  proceeds  from  the  point  of  view  that  the  antenna  is 
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C .  {Continued) 


a  boundary-value  problem  which  can  be  formulated  in  general  terms  and  the  resulting 
equations  solved  by  the  use  of  the  numerical  methods  and  high-speed  digital  computers. 

Thus,  to  solve  this  problem,  equations  are  needed  in  terms  of  parameters  such  as  the  radius 
and  length  that  characterize  the  dipole  antenna,  rather  than  in  terms  of  characteristic  imped¬ 
ance,  terminal  impedance,  etc. ,  which  arc  basically  foreign  to  the  antenna. 

The  dipole  antenna  has  been  investigated  from  the  point  of  view  of  a  boundary  value 

1  7  5  8  9 

problem  by  many  workers  such  as  R.W.  P.  King,  Fallen,  L.  V.  King,  Rrillouin,  and 
Aharoni.1**  In  nearly  all  of  these  papers,  different  but  comparable  methods  were  used. 

Some  are  analytically  complicated  in  determining  ine  equations,  and  ail  are  complicated  in 
their  attempts  to  solve  the  resulting  integral  aquation  Actually,  the  problem  can  be  set  up 
formally,11  and  this  will  lead  to  the  integral  equation  first  obtained  by  Fallen  as  shown 
below. 

1.  The  Differential  Equation 

Since  we  are  concerned  with  antennas  that  will  be  used  at  high  frequencies,  we  will 

assume  the  antennas  to  be  thin.  By  this  it  is  meant  that  the  ratio  of  radius  r  of  th  cross 

o 

section  to  length  21  of  the  antenna  is  negligible.  It  is  not  permissible  to  regard  rQ  as 
zero,  for  the  quantity 


cannot  be  neglected.  Thus,  while  errors  of  the  order  r  /-t,  can  be  tolerated,  the  antenna 

°  9 

cannot  be  regarded  as  infinitely  thin .  A  lower  bound  on  R  has  been  given  by  Brillouin  of 
R  >  14  or  -t-/ro  >  1000  .  For  the  purposes  of  this  presentation  this  is  not  a  significant  re¬ 
striction.  This  restriction  allows  us  to  ignore  the  ends  of  the  antenna,  which  greatly  sim¬ 
plifies  the  problem  but  does  not  decrease  the  useful  accuracy. 

The  analytical  problem  of  determining  the  distribution  of  current  on  an  antenna  which 
we  will  assume  to  be  cylindrical  (triangular  or  other  shapes  can  be  represented  by  a  suit¬ 
ably  selected  cylinder),  of  half  length  l  ,  and  radius  rQ  ,  may  be  formulated  in  terms  of 
the  general  boundary  condition  which  requires  continuity  of  the  tangential  component  of  the 
electric  field  across  any  boundary  surface  between  two  media.  For  this  analysis  it  is  con¬ 
venient  to  use  a  cylindrical  coordinate  system  with  the  z  axis  lying  along  the-  axis  of  the 


C.l  (Continued) 

dipole,  and  z  r  0  at  the  center  of  the  dipole.  We  thus  have  (assuming  the  antenna  to  be  a 
solid  wire  for  purposes  of  discussion)  the  following  boundary  conditions. 

a.  E  as  one  approaches  the  surface  of  the  wire  from  inside  the  wire  must  be  equal 

to  Ez  as  one  approaches  the  surface  of  the  wire  at  the  same  point  from  outside 

the  wire. 

b.  E  as  one  approaches  the  end  of  the  wire  from  inside  the  wire  must  be  equal 

to  Ep  as  one  approaches  the  same  end  of  the  wire  at  the  same  point  from  outside 

the  wire. 

Since  the  ends  of  the  wire  are  required  to  be  small  by  the  restriction  on  SI  as  defined  in 
(a)  above,  the  average  electric  field,  E^  ,  on  the  ends  must  be  less  than  the  average  field 
along  the  wire.  This  follows  since  the  radial  current  must  vanish  at  r  =  0  and  r  =  ro 
and  it  cannot  reach  a  large  amplitude  in  this  extremely  short  distance  relative  to  a  wave¬ 
length.  Accordingly,  nothing  of  significance  is  neglected  as  far  as  the  antenna  as  a  whole 
is  concerned  if  no  account  is  taken  of  the  end  faces  of  the  wire,  and  hence,  Ef  .  We  will 
thus  assume  the  current  to  vanish  at  z  =  ±1  without  flowing  radially  inward  on  the  end 
faces . 

While  at  this  time  conductors  with  finite  conductivity  have  not  been  considered  in  the 
computer  programs,  it  is  planned  in  the  future  to  consider  conductors  with  finite 

conductivity.  Thus,  in  the  formulation  of  the  differential  equation  here,  we  will  assume 
that  the  wire  has  a  finite  internal  impedance  per  unit  length  of  Z.  and  that  the  cross 
sectional  distribution  of  current  in  the  wire  is  essentially  the  same  as  in  an  infinitely  long 
wire.  In  other  words,  for  the  purposes  of  this  analysis  it  will  be  taken  for  granted  that  the 
cross  sectional  and  axial  distributions  of  current  are  mutually  independent.  This  is  always 
true  to  a  very  high  degree  of  approximation  in  a  good  conductor. 

Accordingly  we  have 


E 

z 


r^r 


Z. 

1 


I 

z 


(14) 


where  E  is  the  value  from  the  internal  approach  to  the  surface  and  I  is  the  total  current 
z  z 

in  the  wire  at  point  z  . 
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C.l  (Continued) 

At  any  point  in  space  outside  the  wire  the  electric  field  is 

E  =  -70-jccA  ,  (15) 

where  0  and  A  are  the  scalar  and  vector  potentials,  respectively,  and  co  is  the  angular 
frequency.  Now  the  relationship  between  0  and  A  is 

7-A  =  -jcone  =  j^-^0  •  U6) 

_7 

where  /i  is  the  permeability,  which  for  free  space  has  the  value  4n  x  10  henries  per 

-12 

meter;  e  is  the  permittivity,  which  for  free  space  has  die  value  8. 85  x  10  farads  per 
meter;  and  c  is  the  velocity  of  light. 

By  substituting  Equation  (16)  in  Equation  (15)  we  have 

(;7) 

or 

E  =  uk  +  "'“*)]  -  <18) 


Define 


(19) 


where  X  =  wavelength. 

Since  the  wire  antenna  carries  only  a  z  component  of  current,  the  vector  potential  will 
have  only  a  z  component.  Thus ,  Equation  (18)  will  become ,  with  the  definition  given  in 
Equation  (19), 


(20) 
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If  we  now  impose  the  boundary  condition,  Equations  (14)  and  (20)  will  give  the  following 
differential  equation  in  the  vector  potential: 


or 


Z.I 


1  z 


(21) 


(22) 


From  this  last  equation  we  see  that  the  vector  potential  can  be  represented  by  a  one¬ 
dimensional  wave  equation  which  would  be  homogeneous  for  the  perfect  conductor  (Z.  0). 

It  is  readily  verified  that  the  scalar  potential  can  be  represented  by  a  similar  equation. 

The  differential  equation  given  by  (22)  is  a  nonhomogeneous  equation  and  can  possibly 
lie  rewritten  in  the  form  of  an  ordinary  differential  equation  and  solved  by  numerical  methods 
on  a  igital  computer.  In  contrast  to  initial  value  problems,  which  are  characterized  by 
the  fact  that  the  information  given  concerns  all  the  conditions  at  a  given  point,  a  boundary- 
value  problem  such  as  defined  by  Equation  (22)  is  one  which  gives  the  conditions  at  two  or 
more  distinct  points.  In  this  equation  we  are  seeking  we  must  find  a  relationship  between 
the  vector  potential  and  the  z  coordinate  at  all  points  on  the  wire.  This  type  of  problem 
is  often  called  a  "jury  problem”.  In  contrast,  the  initial  value  problem  is  to  find  the  rela¬ 
tionship  at  another  poinL  given  its  value  at  a  given  point;  i.e. ,  as  the  independent  variable 
is  increased  by  a  small  increment.  This  process  is  sometimes  referred  to  as  a  "marching 
problem”.  The  numerical  solution  of  marching  problems  are  much  simpler  than  those  of 
jury  problems;  in  fact,  all  good  computer  centers  will  have  many  subprograms  to  handle  the 
marching  problem,  while  all  jury  problems  must  be  set  up  and  programmed  separately. 

There  are,  of  course,  other  ways  in  which  one  may  solve  the  boundary  value  differen¬ 
tial  equation  as  shown  in  Equation  (22).  One  approach  that  can  be  used  and  will  be  used 
here  is  to  restate  the  problem  as  an  integral  equation. 


C.2  The  Integral  Equation 


Since  Equation  (22)  is  a  nonhomogeneous  equation,  the  general  solution  must  consist 
of  a  complementary  function  and  a  particular  integral.  The  complementary  function  may  be 
written  in  the  form 

cos  kz  +  CgSinkzj  ,  (23) 

where  Cj  and  are  arbitrary  constants  of  integration.  A  particular  integral  has  been 
obtained  and  is 


(24) 


One  may  verify  that  (24)  does  satisfy  (22)  by  substituting  (24)  in  (22). 

This  equation  is  valid  everywhere  including  the  surface  of  the  antenna,  but  will  not  be 
valid  at  z  =  0  since  the  scalar  potential  will  have  a  discontinuity  due  to  the  driving  voltage 
at  that  point.  This  voltage  source  can  be  inserted  using  King's1  slice  generator.  To  allow 
for  this  additional  condition,  the  differential  equation  must  be  rewritten  as  two  equations 
as  follows: 


z 

=  -J-  cos  kz  +  C2  sin  kz  -  Z.  J  I  (z)  sin  k(z-z)  dz J  for  z  >  0 


(20) 


and 


z 

-i-  cos  kz  +  C4  sin  kz  -  Z.  J  I  (z)  sin  k  (z-z)  dz  J  for  z  <  0 


(27) 


Now  the  boundary  conditions  will  determine  the  constants  Cj  ,  ,  and  . 
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Now  from  Equations  (15)  and  (16)  at  r  =  rQ  the  scalar  potential  is 


r=ro 


(28) 


Inspection  of  Equations  (26)  and  (27)  in  view  of  (28)  indicates  that  if  A  is  an  even  function, 

z 

<t>  must  be  an  odd  function,  and  vice  versa.  If  we  let  that  thickness  of  the  slice  generator 
z 

be  26  ,  ihen 


(29) 


and 


u  -  i  <3»> 

where  V  is  the  applied  voltage,  then  the  function  is  odd.  Thus,  A  must  be  even. 

t*  z 

These  conditions  are  represented  in  equation  form  by 

0=0  (31) 

z  -z 


and 


A 

z 


(32) 


In  order  that  (26)  and  (27)  may  satisfy  (31)  and  (32)  for  all  values,  the  following  relation 
must  hold: 


(33) 
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We  can  then  combine  Equations  (26)  and  (27)  into  the  following  form: 


z 

Az  -  -  ^  [ciCoskz  +  C2  sink|z|  -Z.J  I  (z)  sin  k  (z-z)  dz  j 


It  is  readily  shown  that  (34)  is  unchanged  if  -z  is  everywhere  written  for  z 
gral  one  must  replace  the  variable  z'  by  -z'  . 

The  boundary  condition  on  the  scalar  potential  is 


V  '  6-0  [•,-*-,] 


as  we  allow  the  slice  generator  tl  ickness  to  approach  zero.  From  Equation 
write 


dA 

z 

dz 


z=6 


JU) 

~~2  h 


and 


_^z 

dz 


-JC0 


z=-6 


2  v-6  ’ 


so  that 


0i.  -  -0 


.  2  oA 
3  c  _ z 


-6  oo  dz 


z=-6 


and 


V  -  2 


lim  0  _  2  jc^  lim  ^z 

6-0  6  (jj  6—  0  Qa 


z=6 


(34) 

.  In  the  inte- 

(35) 

(33)  we  may 

(3(5) 

(37) 

(38) 

(39) 
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Upon  differentiating  Equation  (34)  with  respect  to  z  and  setting  z  equal  to  6  and  allowing 
6  to  approach  zero,  we  have 


lim  *AZ 
6-0  ?>z 


z^6 


-*C2 

c 


(40) 


so  that  with  k  -  oo/c  we  then  have 


C2  =TV 


(41) 


The  equation  for  the  vector  potential  then  may  be  written 


I  a.  V 

cos  kz  +  y 


sin  k  z 


-  Z. 

l 


( 


J 

o 


I  (z1)  sin  k  (z-z')  dz' 


(42) 


Figuratively  speaking,  let  us  for  the  moment  sit  back  and  see  what  Equation  (42)  is 
saying,  and  what  it  is  that  we  desire.  First,  it  is  an  equation  for  the  vector  potential  on  or 
outside  the  dipole.  Secondly,  it  has  two  unknown  constants,  and  Z.  .  Z.  is  no  concern 
at  the  moment  since  we  know  it  is  the  internal  impedance  per  unit  length  of  wire  that  can  be 
determined  when  we  so  desire. 

Cj  is  an  integration  constant  and  therefore  related  to  a  boundary  condition.  The  only 
boundary  condition  that  has  not  been  used  is  the  condition  at  the  tip  of  the  dipole.  We  know 
that  the  current  must  be  zero  at  the  tip.  That  is,  C..  can  be  determined  by  the  equation 


z=  l 


(43) 


But  this  is  an  equation  for  the  current  while  Equation  (42)  includes  in  the  equation  for 
the  vector  potential.  Therefore,  it  appears  we  must  obtain  a  relationship  between  the  vector 
potential  and  the  current.  Then,  when  we  have  the  equation  relating  the  vector  potential  to 
the  current  and  substituting  Equation  (42),  we  can  determine  by  the  use  of  Equation  (43). 
To  obtain  the  relationship  between  the  current  and  the  vector  potential,  we  must  use  the 
Helmholtz  equation  which  relates  the  vector  potent’ a  1  and  the  current. 
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i 


The  amplitude  of  the  vector  potential  defined  by  Equations  (15)  and  (16)  must  satisfy 
the  Helmholtz  equation 

V*A  +  k2  A  =  Hi  ,  (44) 

where  T  is  the  volume  density  of  current  flowing  in  the  antenna.  The  Helmholtz  integral, 
by  substituting  in  Equation  (44),  can  be  shown  to  satisfy  this  equation.  This  integral  is 


v 


where  v  is  the  volume  of  the  antenna,  dv'  is  an  element  of  v,  i'j  is  the  current  density 
at  dv' ,  and  r^  is  the  distance  between  the  point  whose  cylindrical  coordinates  are 
(r,  e,  z)  where  A  is  to  be  computed  and  the  element  dv'  is  at  (r1  ,  6',  z') .  There 
are  no  radial  components  to  A  ,  since  we  have  assumed  no  radial  components  of  current 
flow.  The  z  component  is 


A 


z 


It  has  been  shown  by  King  and  Harrison 


11 


that  A  . 
z 


evaluated  from 


A 

z 


=  JL 

4n 


dz' 


(46) 


(47) 


where 


R  =  ^(  z-z')2  +  r2 

on  the  surface  of  the  antenna  differs  by  a  negligible  amount  from  A  computed  by  the 

Z 

exact  equation  given  by  (44)  and  (45) .  Here  R  is  the  distance  from  the  point  (r  ,  0 ,  z) 
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outside  the  conductor  where  A  is  to  be  calculated  to  the  center  of  the  element  dz'  at 

z 

z1  on  the  axis. 

Finally,  by  using  Equation  (47)  in  Equation  (42)  we  have 


dz' 


cos  kz  -*  y  sin  k  |z 


z 

W-*i  /*i 


sin  k  (z-z')  dz' 


(48) 


The  solution  of  this  equation  or  the  same  equation  where  Z.  has  been  set  to  zero  can  be 
done  by  at  least  two  methods .  The  method  we  will  not  use  consists  of  expanding  the  integral 
on  the  left  into  two  integrals  by  the  use  of  integration  by  parts.  This  results  in  one  integral 
which  is  integrable  in  closed  form  and  another  integral  which  is  not  integrable  in  closed 
form.  If  one  performs  this  integration  and  substitutes  back  in  the  original  equation,  a  first 
order  solution  is  derived.  If  Equation  (43)  is  then  applied  and  C1  determined,  one  has 

t  ■*“ 

the  expression  that  was  originally  derived  by  Hallen.  His  derivation,  of  course,  was  not 
the  same  as  given  here.  It  is  this  approach  which  gives  rise  to  the  expression  for  12  as 
used  in  (13).  The  other  method  is  numerical,  and  since  it  is  the  method  used,  the  solution 
will  be  given  in  detail  in  the  next  subsection. 
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D*  NUMERICAL  SOLUTION  OF  THE  INTEGRAL  EQUATION 


At  this  time  we  have  not  considered  conductors  with  finite  conductivity;  therefore,  in  this 
subsection  we  will  assume  that  the  internal  impedance  per  unit  length  is  zero.  We  may  then 
rewrite  Equation  (48)  as  follows: 

e  jj^  dz'  -  B  cos  kz  + 

R 
-l 

where 

4irC9 

B  =  — — — 
jcp 

and 

n  ~  impedance  of  free  space  =  120ir. 

Further, 

e  =  G'  (z,  z'),  (50) 

R 

which  is  the  free -space  Green's  function. 

We  may  then  write  (49)  as  follows: 

l 

U  G'  (z,  z')  dz'  -  B  cos  z  sin  k  j  z  |  =0.  (51) 

-l 

Equation  (51)  is  an  equation  of  the  form 

b 

a(x)y(x)  =  F(x)  +  D  /  K(x,C)  y(4)  d4,  (52) 

a 

where  a,  F,  and  K  are  given  functions  and  D,  a,  and  b  are  constant.  This  equation  is 
known  as  a  Fredholm  integral  equation.  The  function  y(x)  is  to  be  determined. 

When  a  4  0 ,  the  above  equation  involves  the  unknown  function  both  inside  and  outside  the 
integral.  In  the  special  case  when  a  =  0,  the  unknown  function  appears  only  under  the 
integral  sign,  and  the  equation  is  known  as  an  integral  of  the  first  kind;  in  the  case  where 

n-16 


iV 

sin  k  z  =  0 , 


D.  (Continued) 

a  ^  1 ,  the  equation  is  said  to  be  of  the  second  kind.  It  would  appeal*  at  first  glance  that 
Equation  (51)  is  of  the  first  kind,  but  it  is  in  fact  of  the  second  kind.  The  reason  for  this 
is  that  the  integration  constant  B  is  a  function  of  the  current  -  specifically,  the  current  at 
the  tip  of  the  dipole . 

I9 

It  can  be  shown  “  that  a  Fvedhom  integral  equation  can  be  considered  as  the  limit  of  a 
set  of  n  algebraic  equations,  as  the  number  of  equations  increases  without  limit.  Use  can 
be  made  of  this  fact  to  obtain  approximate  solutions  to  such  integral  equations.  Therefore, 
we  must  determine  a  method  whereby  Equation  (52)  can  be  expressed  as  a  set  of  n  algebraic 
equations . 

In  elementary  integral  calculus  an  integral  of  the  form 
b 

y  J  F(x)  dx  (53) 

a 

is  defined  as  the  limit  of  an  equation  of  die  form 


y 


lim 
n  -® 


*Xk  F<Xk>’ 


(54) 


where  the  interval  (a,b)  is  divided  into  n  subintervals  of  lengths  ux. ,  Ax„ . Ax  ,  and 

jl  u  n 

Xj.  is  the  point  in  the  k-th  subinterval.  An  approximate  value  of  y  can  be  obtained  by  not 
proceeding  to  the  limit,  and  hence  b>  expressing  y  approximately  as  the  weighted  sum  of 

,  x2 . xr  of  the  interval  (a,b): 


the  ordinates  F(x^ )  at  n chosen  points  x^ 


n 


h  k?,  wk  F<\> 


(55) 


where  h  is  the  length  of  all  intervals  and  Wj,  is  a  "weighting”  coefficient  associated  with  the 
point.  In  its  simplest  form,  we  let  the  point  x^  be  at  the  center  of  the  k-th  interval 
and  give  the  value  unity  to  all  the  weighting  coefficients.  The  resulting  formula  is  known  as 
the  trapezoidal  rule.  The  accuracy  of  the  approximation  can  of  course  be  increased  by 
increasing  n  or  changing  the  weighting  coefficient  or  functions.  An  example  of  the  latter 
is  Sunpson's  rule,  which  is  applicable  only  if  n  is  odd,  is  where  the  weighting  coefficients 
are  of  the  form 


1 1  =  1/3 jl. 


4,  2,  4 . 2, 
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13  14  1 5 

when  n  is  greater  than  five.  More  elaborate  formulas  can  be  found  in  the  literature  ’ 
where  this  technique  is  often  called  Gaussian  quadrature. 

Thus,  Equation  (52)  by  the  use  of  Equation  (55)  can  be  approximated  in  the  form  when  it 
is  the  second  kind  where  D  =  1 .  Let  y(x)  not  under  the  integral  be  B(x): 


B(x)  +  F(x)  +  h 


n 

E  W.K(x,xJ  y(x.) 
k=l  ‘v  K  K 


0 


(57) 


where  the  points  are  spaced  a  distance  h  apart  with  the  first  and  last  points  being  h/2 
from  the  points  a  and  b.  If  we  now  require  that  equation  be  written  at  each  of  n  chosen 
points  where  the  points  are  the  x^  points ,  we  obtain  n  linear  equations 


n 

B(x.)  +  F(x.)  +  h  £ 
k=l 


WkK(xi,Xk)  y(xk) 


0 

i  =  1,  2, . . .  .n 


(58) 


in  the  unknowns  y(x  ),  y(x„), . . .  .y(x  )  which  specify  approximate  values  of  the  unknown 
I  a  n 

function  y(x)  at  the  n  points.  Introducing  the  following  abbreviations , 


F. 

l 


=  F(x.),  K..  =  K(xr  x.), 


B:  =  B(x.) , 

X  1 


where  K„  is  the  value  of  K(x, £)  when  x  =  x^  and  £  =  x^,  the  set  of  equations  given  by  (58) 
can  be  written  in  the  form 


hk?1Kik>'k  +  FrBi  =°- 


(59) 


Thus,  if  we  consider  the  numbers  F.,  B.,  and  y,  as  components  of  the  vectors  F,  B,  and 

IX  K 

y ,  and  define  the  matrix  K  =  [K„] ,  the  set  of  equations  (59)  can  be  written  concisely  in 
the  form 


h  X  [y]  +[FJ  +  [B]  =  0 


(60) 


or 


[y]  =l/h  K 


-1 


[Fl  +  [B  ] 


=  0, 


(61) 
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where  K  *  =  inverse  of  the  matrix  R,  thus  obtaining  the  value  of  y  at  n  uniformly  spaced 
points. 

Using  the  method  outlined  in  the  above  paragraphs,  we  may  rewrite  Equation  (51)  as  a 
set  of  linear  algebraic  equations 


h  X!  wk  G'  <zi>  V  !k  "  B  C0S  k2!i  +irj  sin  k  j  zi  |  *  “  i»2, . .  .n, 
k-1 


where  Ij,  is  the  current  on  the  dipole  at  the  k-ih  point .  Investigation  of  this  set  of  n  euua- 
tions  will  show  that  we  have  n  equations  and  n+2  unknovns;  the  unknowns  being  B  and  V. 
We  must  find  two  additional  equations  before  we  can  solve  the  set  of  equations  given  in  (62) . 
The  applied  voltage  which  is  represented  by  V  presents  no  problem .  We  may  let  the  applied 
voltage  be  unity  giving  the  following  equation: 


V  -  1  =  0. 


(63) 


To  determine  B  we  must  write  an  equation  that  expresses  the  current  at  the  tip.  This 
can  be  accomplished  in  two  ways.  The  first  method  we  used  was  to  write  a  Taylor  series 

expansion  of  the  current  at  the  points  k  “  1,  2,  3, - m  that  is  in  terms  of  1^,  I9,  Ig, - 

I  ,  where  m  <  n,  since  there  are  only  n  points  on  the  half  dipole.  The  k  =  1  point  is  h/2 
from  the  tip,  k  -  2  point  is  3h/2  from  the  tip,  k  =  3  point  is  5h/2  from  the  tip,  and  the 
k  =  m  point  is  (2m  +  l)h/2  from  the  tip.  While  the  resulting  equation  does  not  involve  B, 
it  is  the  additional  expression  required  to  give  the  required  square  matrix.  Since  we  are  no 
longer  using  the  Taylor  series  expression,  the  lengthy  derivation  will  not  be  given,  only 
the  final  equation  which  is  as  follows: 


where 


is  the  binominal  coefficient  given  by 


/m\_  ml 

V  3/  )■  (n  "  3)'. 


(65) 
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The  second  method  is  to  extrapolate  the  currents  1^ ,  I2 ,  and  1^  which  are  unknowns  and 
require  that  the  current  at  the  tip  be  zero .  The  equation  will  be  of  the  form 

f  (l  -  h/2) Ij  +  g(t  -  I2  +h  I3  =  0.  (66) 

The  form  of  f,  g,  and  h  will  be  determined  later.  If  the  set  of  equations  represented  by 
Equation  (62)  are  applied  over  both  halves  of  the  dipole ,  as  is  inferred ,  there  are  two  things 
in  error.  First,  with  numerical  integration  represented  by  (62)  as  with  analytical  integra¬ 
tion,  one  cannot  integrate  over  a  discontinuity.  As  shown  in  Equation  (35),  there  is  a  dis¬ 
continuity  in  the  scalar  potential  at  the  slice  generator  in  the  center  of  the  dipole.  Secondly, 
since  the  current  distribution  on  the  dipole  is  an  even  function  and  Equation  (62)  is  identical 
for  z  and  ,  the  resulting  matrix  will  be  singular.  These  two  things  can  be  corrected 
quite  simply,  while  at  the  same  time  decreasing  the  number  of  terms  in  the  matrix  (62)  by 
two,  by  taking  note  of  the  symmetry  of  the  Green’s  function  which  is 

G(z.,  z.)  =  G(z.,  z.).  (67) 

First  we  reduce  the  number  of  points  from  2n  to  n:  that  is,  we  integrate  only  over  one  half 
of  the  dipole.  At  the  same  time  we  replace  the  free-space  Green's  function  with  a  modified 
Green's  function  which  includes  the  effects  of  the  point  -z..  That  is,  let 

G  *zi’  zk*  =  G'  *zi’  zk*  +  G'  *~zi’  zkJ'  (68) 

Thus,  by  writing  the  set  of  equations  represented  by  Equation  (62)  using  the  modification 
given  by  Equation  (68) ,  Equation  (63) ,  and  Equation  (64) ,  or  Equation  (56) ,  we  have  a  set  of 
algebraic  equations  that  relates  the  current  on  a  dipole  to  the  length  of  the  dipole,  the  radius 
of  the  dipole,  and  the  applied  voltage. 

E.  SINUSOIDAL  APPROXIMATION  FOR  THE  NUMERICAL  SOLUTION 

In  subsection  D  it  was  indicated  that  more  elaborate  formulas  for  the  weighting  coeffi¬ 
cients  can  be  found  in  the  literature.  Most  of  these  methods  are  written  assuming  one  has 
no  knowledge  of  the  behavior  of  the  unknown  function  under  the  integral  sign.  Dr.  Mei 
has  pointed  out  that  we  do  have  knowledge  of  the  behavior  of  the  current  on  the  dipole .  We 
know  that  we  must  be  able  to  represent  it  by  a  Fourier  Series  of  finite  length.  Further, 
we  should  be  able  to  represent  the  current  in  a  region  ±h/2  about  a  point  quite  accurately 
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E .  (Continued) 

by  only  the  first  three  terms  in  a  Fourier  series.  Thus,  we  may  approximate  the  current 
i (■/.)  near  the  k-th  segment  as 


I(z)  =  +  Bk  sin  kfz-z^)  +  cos  k(z-zk),  (69) 

or  by  applying  at  the  k-th  point  and  the  points  on  either  side  and  writing  Ik  for  I(zk),  we 


have 

\-l 

=  sin  kh  +  Ck  cos  kh 

(70) 

*k 

*  \  +ck 

(71) 

Vi 

=  A^  +  B^  sin  kh  +  C-k  cos  kh 

(72) 

Solving  this  set  of  simultaneous  equations  for  A^,  B^,  in  terms  of 

W  V  W,  wchave 

=  2  sin  kh  cos  kh  -2  sin  kh  -  2  sin  kh  (cos  kh-1) 

(73) 

Ak 

=  |-sin  kh  +  2  sin  kh  cos  kh  1^  -  sin  kh  I^J 

(74) 

Bk 

=  ~  |(l-cos  kh)  Ik  l  +  (cos  kh  -  cos  kh)  ^  +  (cos  kh-1)  Ik+1  j 

(75) 

Ck 

‘  k  |sin  kh  ‘k-l  ‘2  sin  kh  ^  +  sm  kh  ^,1  . 

(76) 

This  gives  for  I(z) 

l(z) 

=  ~  |  Ik  l  |-sin  kh  +  (1-cos  kh)  sin  k(z-zk)  +  cos  k(z-zk)  sin  kh  J 
+  Ik  (2  sin  kh  +  cos  kh  -  2  sin  kh  cos  k(z-zk)  j 
+  Ik+1  |  -sin  kh  +  (cos  kh-1)  sin  k(z-zR)  -  sin  kh  cos  k  (z  -  zk)J  |  . 

(77) 

E .  (Continued) 

Writing  Equation  (77)  in  a  more  compact  form 
I(z)  =  X(z)  IR_X  +  Y(z)  ^  +  Z(z)  Ik+1 


where  the  coefficients  are  given  by 


X(z)  -  ^  J  -sin  kh  -*  (1-cos  kh)  sin  k  (z-zk)  +  sin  kh  cos  k(z-zk)  J  (79) 

Y(z)  =  ^  1 2  sin  kh  cos  kh  -2  sin  kh  cos  k  (z-zk)  j  (80) 

Z(z)  =  J-sin  kh  +  (cos  kh-1)  sin  k  (z-zk)  +  sin  kh  cos  k  (z-zk)|  .  (81) 

Equation  (78)  states  that  the  current  in  a  given  segment  is  a  function  of  the  current  at  the 
center  of  that  segment  and  the  current  at  the  center  of  the  segments  on  either  side .  Thus , 
as  it  stands  (78)  cannot  be  used  for  the  point  k  =  1  located  h/2  from  the  tip,  or  the  point 
k  =  n  located  h/2  from  the  center,  since  these  points  do  not  have  points  on  either  side. 
Therefore,  we  must  determine  a  special  set  of  coefficients  for  each  of  these  points  using 
the  next  two  points  moving  away  from  the  boundaries.  Consider  the  tip  point  first.  The 
expression  for  I(z)  in  the  tip  segment  is: 

I(z)  =  +  Bj  sin  k  (z-Zj)  +  cos  k  (z-Zj)  (82) 


X1  =  A1 


1^  =  sin  kh  +  cos  kh 


Ig  =  Aj  +  Bj  sin  2  kh  +  cos  2  kh. 


Solving  this  set  of  simultaneous  equations  we  have 


*1  =*k 


Aj  =  —  j -sin  kh  Ij  +  sin  2  kh  I2  -  sin  kh  Ig 
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B1  "  ^  [<cos  **  -  cos  2kh)  Ix  +  (cos  2kh-l)  Ig  +  (1-cos  kh)  Ig] 
^1  -  [  (s“v  ^kh  -  sin  kh)  1^  -  sin  2kh  Ig  +  sin  kh  Ig  j  . 


Thus, 


«Z>  -  x1i1+y1i2  +  zii3 


for  the  tip  segment  where 

X1  =  [-sin  kh  +  (cos  kh  -  cos  2kh)  sin  kfz-Zj) 

+  (sin  2kh  -  sin  kh)  cos(z-z]L)  | 

[sin  2kh  -  (1-cos  2kh)  sin  k(z-z^)  -  sin  2  kh  cos  k(z-z^)  J 

Z1  =  ^  [ _sin  ^  +  (1-  cos  kh)  sin  k(z-Zj)  +  sin  kh  cos  k(z-zx)  j 
The  expression  for  the  n-th  segment  is 

I(z)  =  An  +  Bn  sin  k(z-zn)  +  Cn  cos  k(z-zn)  , 

which  gives  using  the  above  method 


where 


I(z)  =  X  I  +  Y  I  +ZI 

n  n  n  n-1  n  n-2’ 


xn  ::  ^  [-sin  kh  +  (cos  2kh  -  cos  kh)  sin  k(z-z  ) 
n  n 


(sin  2kh  -  sin  kh)  cos  k(z-zn) 


(88) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 
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Yr  =  (sin  2kh  +  (1  -  cos  2kh)  sin  k(z-z  )  -  sin  2kh  f>os  k(z-z  )j  (97) 

n  n  ' 

Z  =  J  -sin  kh  +  (cos  kh-1)  sin  k(z-z  )  +  sin  kh  cos  k(z-z  )j  .  (98) 

n  n  * 

In  Equation  (55)  it  is  assumed  that  the  function  was  constant  over  the  segment,  thus  the 
area  was  represented  by  the  product  of  f(xk)  and  h.  Now,  however,  I(z)  is  not  constant  and 
to  obtain  the  area  we  must  perform  an  integration  over  the  distance  h.  Further,  each 
element  of  that  part  of  the  matrix  which  represents  the  integral  will  consist  of  terms 
involving  the  currents  on  either  side  excluding  the  special  cases  near  the  boundaries. 

These  terms  of  the  matrix  will  be 

bk 

A(i,k)  =  f  (XkIk_1+YkIk  +  ZkWG(z.,z-)dz-,  (99) 

\ 

where 

ak '  \  - h/2 

bk  *  zk  +  h/2 
and  for  the  tip  point 


A(i,  1)  =  J  (Xj  Ix  +  Y  I2  +  Z  I3)  G  (z.,  z«)  dz'  ,  (100) 

al 

where 

ax  =  zx  -  h/2  =  l  --Y 

bj  =  zx  +  h/2  =  l, 
and  for  the  base  point 

b 

n 

A(‘  •„)=/  <\I„tynIn-l+ZnWG<VZ'>dZ’  |I01) 
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where 


a  =  z  -  h/2  =  0 
n  n 


b 


n 


z  +•  h/2 

n 


3h 
2  ' 


Since  the  unknown  currents  1^ 
ask  if  replacing  one  integral  equation  by  m  integral  equations  do  we  represent  the  integrals 

'  V  and  w 

are  not  variables  in  the  integrals  -  they  are  constants.  The  type  of  integral  we  are  con¬ 
cerned  with  here  is  of  the  form 


in  (82)  to  (84)  by  a  m  point  sum.  This  is  not  true  since  the  current  L 


-1’ 


1^,  and  are  under  the  integral  sign,  one  might 


b 

Jx  (z\  zR)  G  (z.,  z’)dz*  (102) 

a 


where  b-a  =  h.  When  given  z',  and  z.  we  can  calculate  the  product  XG 

One  of  the  more  accurate  methods  of  numerical  integration  is  using  Lagrangian 
16 

integration  coefficients ,  °  where 
b 

J  f(x)  dx  =  aQ  yQ  +  y1  + - <*n  yR,  (103) 

a 

where  the  ot's  are  the  coefficients  and  y's  are  the  values  of  f(x)  at  the  n  points.  Accuracy 
increases  with  the  number  of  points,  but  the  computer  time  also  increases.  The  use  of 
10  intervals  which  is  11  points  seems  to  be  a  good  compromise.  The  value  of  the 
Lagrangian  integration  coefficients  for  n  =  10  are  given  in  Table  III— 1 .  Thus,  by  evaluating 
the  product  X.G.a.  at  each  of  the  i  points  in  the  interval  (a,b)  and  summing  and  multiplying 
by  h/10,  the  integral  represented  by  Equation  (102)  can  be  evaluated  and  inserted  in 
Equation  (99) -(101). 
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Table  1II-1 .  Lagrangian  Integration  Coefficients 


i 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0.268341483 

1.775359414 

-0.810435705 

4.549462882 

-4.351551226 

7.137646304 

-4.351551226 

4.549462882 

-0.810435705 

1.775359414 

0.2683414o3 


E .  (Continued) 

With  this  sinusoidal  approximation  we  can  decrease  the  number  of  points  on  the  dipole 
for  a  given  accuracy.  Also  we  need  no  longer  do  a  Taylor  series  expansion  for  the  boundary 
condition  but  can  use  Equation  (90)  with  z  =  l. 

I  (t)  =  Xx  (l)  Ij  +  Yx  (l)  I2  +  Zx  (l)  I3  =  0  (104) 

and  the  current  at  the  feed  point  is  given  by 


I  (0)  =  X  (0)  I  +  Y  (0)  I  +  Z  (0)  I 
'  '  n  '  ’  n  nv  '  n-1  n'  •  n- 


(105) 


F.  SIMULTANEOUS  LINEAR  EQUATIONS 

Finally,  one  note  on  the  solution  of  the  n  +  2  simultaneous  linear  equations.  Most 

computing  centers  have  subroutines  that  can  solve  a  set  of  complex  simultaneous  linear- 

equations  using  various  methods.  While  the  computer  we  have  available  is  sufficiently 

large  to  handle  a  single  dipole,  it  was  felt  that  when  one  had  antennas  with  more  than  five 

dipoles  there  would  be  insufficient  core.  Thus,  a  subroutine  was  developed  that  does  not 

require  more  than  two  rows  of  the  argumented  matrix  in  the  core  at  one  time  by  the  use  of 

magnetic  tape  units  for  storage.  The  method  used  was  the  Gauss-Jordan  elimination 

17 

method  with  normalization.  This  method  allows  us,  if  desired,  to  solve  up  to  500  simul¬ 
taneous  linear  equations  without  overflow.  By  properly  filling  the  matrix  we  have  not 
experienced  any  difficulty-  with  round-off  errors  so  far  (130  equations). 
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G.  EQUATION  FORMULATION  FOR  A  LOG-PERIODIC  DIPOLE  ANTENNA 


After  one  has  formulated  the  equations  for  a  single  dipole,  what  changes  or  additions 
must  be  made  to  give  the  set  of  simultaneous  linear  equations  that  will  describe  a  log- 
periodic  antenna?  There  is  one  change  to  be  made  and  one  addition  to  be  made.  We  must 
change  Equation  (t32)  so  that  the  interaction  between  the  different  dipoles  of  the  antenna  are 
properly  handled.  We  must  also  add  additional  equations  to  represent  the  transmission 
lines  that  tie  the  dipoles  of  log -periodic  elements  together. 

The  approach  used  below  to  determine  the  changes  is  based  on  the  separate  work  of 
one  of  the  authors .  ^ 

The  vector  potential  on  a  dipole  antenna  in  the  presence  of  N-l  other  dipoles,  which 
need  not  be  of  the  same  length,  can  be  represented  by  an  equation  which  is  based  on 
Equation  (47)  using  the  symbology  given  in  Equation  (50).  The  vector  potential  will,  if  all 
the  dipoles  are  parallel  to  the  z-axis  with  their  centers  being  in  the  z  =  0  plane,  have  only 
a  z  component.  Thus,  we  have  at  point  P 

,  1 

nd  --m 

Az(p)  "  E  /  I  (z'>  G  (z\  z^  dz*  ,  (10G) 

m=l  ~ 


where  A  (P)  is  the  vector  potential  at  any  point  which  is  outside  the  dipoles,  nd  is  the 
z 

number  of  dipoles,  is  the  length  of  one  half  of  dipole  m,  Im(z')  is  the  current  flowing 
in  dipole  m  at  the  point  z\  and  G(z\  z^)  is  a  modified  Green's  function  of  the  form  given 
by  Equation  (68)  and  with  r  of  Equation  (47)  replaced  by  d,  the  horizontal  distance  between 
the  dipole  m  and  point  P . 

The  scalar  potential  at  point  P  can  be  found  from  the  following  relationship 
-jkR 


0  = 


/ 


oe 


4ff«R 


dv, 


(107) 


where  p  is  the  charge  density  and  v  is  the  volume.  We  may  rewrite  Equation  (107)  for  our 
antenna  as  follows: 


0<P) 


d*m  l  G  (z1,  zp)  dz'  , 

az' 


(108) 
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which  after  integration  by  part  and  using  the  quantity 


we  have 


nd 


__I_  V  /  I  fzM  — L  G<2’.  zp)  dz' 
jwC  n£l  J  m  &P 


The  tangential  electric  field  on  any  antenna  produced  by  the  currents  on  the  other  antennas 


Ez  -  -  -  jwp  Az  =  ^P)  -  jwp  A^fP), 


where  the  point  P  is  on  the  surface  of  the  dipole  but 

,  3A  (P) 

__  ^ 


jwc  oz. 


Therefore, 


1  3“Az(P) 


•1WEZ  =  SI"  +  k  Az.<P>  • 


For  the  tangential  electric  field  to  vanish  on  the  surface  of  the  dipoles,  we  have  for  dipole  m 

E  (P)  -  -E  (P)  (incident)  =  -V  6(1  )  ()  io) 

zv  'm  z'  'm  '  m  '  m  ' 


A  C,  sin  Viz  +  C„  cos  k  z 
zp  1  1  p  2  p 
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This  last  equation  can  now  be  represented  by  a  set  of  algebraic  equations  in  the  form  given 
by  Equation  (62) ,  etc . , 

nd  l  n  j  2V. 

£  \  h  51  W  I.  G(z.  ,  z.)  >  -  3.  cos  kz.  +  T"T^  sink  z.  =  0  (117) 

.  ,  )  ,  m  lm  '  un  j  t  j  i  27)  t  '  ' 

1=1  l  m-1  J  J  J  o 

j  =•  1,  2, - n. 

Here  we  now  have  a  boundary  condition  for  each  dipole,  which  is  the  reason  that  the  B's  are 
subscripted,  and  each  dipolj  has  its  own  driving  voltage.  Also,  while  Equation  (117)  cannot 
show,  as  it  is  written,  the  use  of  the  sinusoidal  approximation  is  used  as  outlined  in  sub¬ 
section  E.  The  equations  for  the  boundary  conditions  on  the  tips  of  the  dipoles  are  written 
using  the  form  given  by  Equation  (104) . 

H.  TRANSMISSION  LINE  EQUATIONS 

Equation  (63)  which  is  the  single  equation  for  the  driving  voltage  on  a  single  dipole  will 
become  the  driving  voltage  of  the  transmission  line  which  feeds  the  log-periodic  dipole 
elements.  To  relate  this  voltage  to  the  base  of  each  dipole,  an  additional  set  of  equations 
in  terms  of  the  current  at  the  base  of  each  dipole,  flowing  into  the  dipole  and  down  the  trans¬ 
mission  line  must  be  determined .  These  equations  are  based  on  the  assumption  that  the 
transmission  line  has  no  ohmic  or  radiation  losses. 

There  is  no  need  to  derive  or.e  single  equation  that  gives  the  relationship  between  the 
injait  voltage  or  current  and  the  voltage  or  current  at  the  most  distant  dipole.  We  need  only 
derive  the  relationship  between  the  n  and  n+l  dipole  and  repeat  this  relationship  between 
the  nd  and  the  n+2  d’poles  until  all  sections  of  the  transmission  line  have  been  covered, 
rigute  III-l  shows  the  geometry  involved. 

19 

The  equations  for  an  ideal  transmission  line  are  given  by 

V  -  Ve"-*sn  +  V'ejsn  (118) 

n 

Z  i  -  Ve-jsn  -  V’cJsn  (119) 

on 

where  is  the  distance  to  a  reference  point  on  the  left  in  radians  from  point  n,  V  is  the 
amplitude  of  the  direct  wave,  V'  is  the  amplitude  of  the  reflected  wave,  V  is  the  voltage 


Figure  II- i.  Transmission  Line  Geometry 
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across  element  n,  in  is  the  current  flowing  into  the  junction  of  the  transmission  line  and 
the  r>  element,  while  ZQ  is  the  surge  impedance  of  the  line.  The  determinant  of  these 
two  equations  is 


D 


-js  js 
e  J  n  eJ  n 


e  sn  -ejSn 


=  -2 


while 


e_jSn  V 


n 


-js  Z  i 

T 


...  i  e  J  n  on 
v  - - 


“JS 
e  J  n 


(V  -  Z  i  ) 
'  n  on' 


and 


V  = 


V  eJSn 


„  .  js 
Z  1  -eJ  n 
o  n 


js 

-e  n 


(V  +  Z  i  ) . 

'  n  r\  tv 


-2  n  on 


(120) 


(121) 


(122) 


Now  let  us  consider  the  equations  at  the  n+1  dipale  with  reference  to  the  n  dipole.  If 
the  spacing  between  the  dipoles  were  sufficiently  short  so  that  the  effects  of  the  transmission 
line  were  zero,  we  need  only  consi  er  the  transposition  of  the  wires.  Thus, 


V  ..  «  -V  or  -V  ,  i  =  V 
n+1  n  n+1  n 


(123) 


Z  i  =  -Z  i  or  -Z  i  . 
o  n*l  on  o  rr 


Z  i 
o  n 


(124) 


or  rewriting  (118)  and  (119)  fcr  the  n+1  dipole  we  have 


-V  -  Ve"JSn+l  +V'ejSn+l 
n  U 


(125) 


-Z  i'  -  Ve  J“n+1  -VV“n+l  . 
o  n-*l 


(120) 
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Substituting  Equations  (121)  and  (122)  in  (125),  we  have  after  reduction 


-V  =  v  |  ejP-n  *  e  jPn  1  +  JZ  i  e^n-  ~  e~^n 

-  1  2  J  on  2j 


n+1  n 


(127) 


where 


( T  =  S  -  S  , 

n  n  n+1. 


which  gives 


-V  . ,  =  V  cos  a  +  jZ  i  sin  a 
n-^1  n  n  J  o  n  n 


(128) 


or 


V  ,  V  cos  o  +  jZ  i  sin  a  =  0 
n+1  r»  n  J  o  n  n 


(129) 


In  a  similar  manner,  we  have 


jvn 

‘U  +  T"  sin(rn  +  j‘noosar  *  0  ' 
o 


(130) 


The  expression  for  the  current  i^+1  is  for  the  current  to  the  right  of  dipole  n+1 ,  while 
the  equation  needed  is  for  the  current  just  to  the  left  of  dipole  n+1 .  The  voltage  across  the 
transmission  line  as  one  moves  toward  the  generator  and  passes  the  feed  point  of  the  n+1 
dipole  must  be  continuous  so  Equation  (129)  is  correct  as  it  stands.  However,  the  current 
must  increase  by  the  amount  flowing  into  the  dipole.  Therefore,  we  must  rewrite  Equation 
(130)  to  include  the  current  I  the  current  flowing  into  the  n+1  dipole.  The  current  on 
the  left  must  be 


*L  ^n+l^L  +  *0+1’ 


while  the  current  on  trie  right  is 


lR  ~  ^+1  ' 


(131) 


(132) 
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Therefore,  we  have 


jV 


i  . ,  -  I  , ,  +  — —  sin  a  ■>  i  cos  <7  =  0  , 

n-rl  n+1  Z  n  n  n 

o 


(1 33) 


where  I  ^  is  determined  by  liquation  (95)  as  applied  to  the  n+1  dipole. 

To  demonstrate  how  we  apply  .hese  equations  and  consider  the  ends  of  the  transmission 
lines,  consider  the  5 -element  dijxlc  shown  in  Figure  III -2.  The  current  flowing  in  an 
antenna  at  the  base  of  the  antenna  is,  using  Equation  (45), 


I.  =-  X.I.,  -!  Y  .,1 .,  ,  +  Z.I.,  o 

j  jjn  j.In-1  kjn-2 


(1*14) 


where  the  single  subscript  on  1  inuicates  the  current  flowing  in  to  the  base  of  the  j  dipole. 
The  equations  for  the  voltages  will  be 

v- 


sin  i 


V,  + 


1  cos  a 


1 


■*  J  i 


1 


l  cos  a i  o 


Z  *  0 


(135) 


V„  sin 

V2  +  cos  o'  ^  *2  cos  a9  Zo  -  0 


(13G) 


V, 


sin  a.. 


V  + 

v  n 


3  cos 


+  3  i. 


3  cos  O' ^  ° 


Z  =  0 


(137) 


Vr  sin  o, 

V ,  +  — +  j  i , - -  Z  =  0 

4  cos  a,  4  cos  o .  o 

4  4 


(138) 


V  sir,  or 

V .  +  — 2™  +  j  i_ - iL  z  -  0 

:>  cos  o.  a  cos  O-  o 
o  5 


(139) 


V  -1.0 
o 


(140) 
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For  the  currents  we  have,  where  i's  represent  transmission  line  currents  and  I's 
represent  dipole  currents, 


‘i  -  £i  =  0 


iVi 

i9  -  I2  +  "2 —  sin  +  cos  a1  =  0 
o 


jV 

i3  ~  *3  +  l1  ^cr2  +  i2cosa2  =  0 
o 


jV3 

i4  ~  !4  +  ~smCr3  +i3COSCT3  =  ° 
o 


jV4 

i_  -  I_  +  -= —  sin  a  +  i.  cos  cr  =  0 
o  a  Z  4  4  4 

o 


jVS 

in  ~  7? —  sin  O-  -  i-  cos  a.  =  0 

0  Z  o  o  o 

o 


(14 1  > 

(142) 

(143) 

(144) 

(143) 

(143) 


It  is,  of  course,  not  realistic  to  have  the  long  transmission  line  in  front  of  the  antenna, 
but  this  has  been  left  on  to  make  it  easier  to  write  the  equations.  It  can  be  removed  by 
subtracting  cr_  from  all  a's. 

I.  FAR  FIELD,  INPUT  IMPEDANCE,  AND  GAIN 

20 

The  relationship  between  the  current  on  a  short  dipole  and  the  far  field  is  given  by 


j  6C v  If,  sin  0 

*9  IX 


(147) 


where  I  is  the  current,  0  is  the  polar  angle,  s  is  the  distance  from  dipole  to  distant  point 
where  Eg  is  to  be  determined,  \  is  the  wavelength  in  meters,  and  l  is  the  length  of  the 
dipole.  It  is  required,  of  course,  that  t«  \  and  that  the  current  be  constant  over  f  or  I 
represents  its  average  value.  The  log -periodic  antenna  as  we  have  represented  it  does  not 
consist  of  a  single  short  dipole,  but  arrays  of  short  dipoles.  Thus,  the  current  shown  in 
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I .  (Continued) 

Equation  (147)  is  only  representative  of  all  the  currents  represented  in  Equation  (117).  We 
need  then  some  relationship  that  will  give  the  total  Eg  in  the  far  field. 

If  we  assume  that  the  distance  to  the  point  P  where  Eg  is  to  be  determined  is  suffi¬ 
ciently  far  that  the  lines  from  P  to  all  of  the  current  elements  on  the  antenna  can  be  con¬ 
sidered  parallel;  that  is,  we  are  in  the  Fraunhofer  region  and  we  can  rewrite  Equation  (147) 
as  follows 

ND  NP  i  6011 1.  sin  8  <Tjk  (sim  '  V 

Ee  “  £.  E  - -  <»» 

m-i  i'l  o 

where  sq  is  the  distance  to  a  reference  point,  for  example,  the  apex  of  the  log  periodic  and 
s.  is  the  distance  to  the  current  element  i.  NP  is  the  number  of  points  on  the  m-th  dipole 
and  ND  is  the  number  of  dipoles  in  the  array 

Let  yim  represent  the  distance  from  the  center  of  the  j  dipole  to  the  i-th  point  in 
radians.  As  before,  @m  represents  the  distance  from  apex  of  the  log -periodic  to  the  m 
dipole  in  radians.  Then,  by  the  use  of  a  little  trigonometry,  we  have 

k  (s.  -  s  )  =  -y.  cos  0+o  sin  0  cos  0,  (149) 

un  o  un  m 

where  0  is  the  polar  angle  and  0  is  the  azimuthal  angle  in  a  spherical  coordinate  system. 

Thus,  the  far  field  E»  at  point  P  at  a  distance  s  ,  polar  angle  0,  and  azimuthal  angle  0, 

o  o 

is  expressed  by 

ND  NP  j  6077  I  sin  0  , 

Efl  =  £  H  - r-> - exp<-j(-y  cos  6  +o  sin  0  cos  0)}  (150) 

0  m=l  i=l  V  (  im  m  f 

If  the  applied  voltage  to  the  antenna  is  one  volt,  then  the  real  part  of  the  inverse  of  the 
input  current  is  the  radiation  resistance  while  the  imaginary  part  is  the  radiation  reactance. 
The  input  current  is  shown  in  Figure  III -2  as  i  .  The  transmission  line  between  the  shortest 
dipole  and  the  apex  is  not  a  part  of  the  antenna,  it  is  assumed  that  the  technique  given  after 
Equation  (146)  is  used  which  makes  this  section  of  line  of  zero  iength. 

Gail,  as  it  will  be  used  on  this  project  is  defined  as  the  gain  relative  to  a  lossless 
isotropic  antenna  in  free  space.  It  will  be  expressed  in  decibels — i.e  — dbi  The 


I.  (Continued) 


mathematical  definition  of  gain  then  is,  which  is  in  agreement  with  the  IEEE  definition  of 
power  gain  in  a  given  direction, 

c  io  '  >  [  x  the  radiation  intensity  in  that  direction  ,  r 

ubi  i0^  I  total  power  accepted  by  the  antenna  ’  ' 

where  the  radiation  intensity,  P,  is  the  power  radiated  in  the  given  direction  per  unit  solid 
angle.  That  is 


P 


2 

r 


120ff 


(152) 


J.  CONCLUSION 

This  completes  the  derivation  and  description  of  the  expressions  that  relate  the  physical 
parameters  of  a  log-periodic  antenna  to  the  currents  on  the  elements  of  the  antenna,  the 
input  impedance,  and  gain  where  the  solution  was  obtained  assuming  that  we  have  a  boundary 
value  prubiem.  Also  the  solution  was  obtained  knowing  that  the  best  numerical  techniques 
and  digital  computers  were  to  be  used. 
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Ill  NUMERICAL  RESULTS 


A.  _ INTRODUCTION 

During  the  study  covered  by  this  report  a  step-by-step  procedure 
was  used  to  obtain  a  computer  program  for  a  log-periodic  antenna  in  free 
space  and  over  real  earth.  At  first  the  theory  presented  in  Section  II  of  this 
report  was  programmed  to  calculate  the  current  distribution  on  just  a  single 
dipole.  When  the  validity  of  the  results  of  this  simple  program  had  been 
verified,  the  more  complicated  case  of  two  Identical  dipoles  in  free  space 
was  mathematically  modeled  and  programmed.  This  procedure  of  applying 
the  theory  to  problems  of  ever-increasing  complexity  and  verifying  the 
calculations  at  each  stage  of  the  process  was  followed  throughout. 

The  following  subsections  present  the  numerical  results  of  some  of 
the  most  interesting  stages  in  the  development  of  the  present  computer  pro¬ 
gram.  At  each  step  the  computer  results  are  compared  with  data  previously 
published  or  to  data  obtained  experimentally  wherever  practical. 

B.  _ SINGLE  DIPOLE 

As  stated  above,  the  theory  was  first  applied  to  a  single  dipole  in 
free  space.  The  dipole  was  divided  into  2N  sections  and  the  current  computed 
for  each  of  the  individual  sections.  At  first  the  current  was  assumed  to  be 
constant  on  each  section  of  the  dipole,  but  later  it  was  found  that  the  number 
of  dipole  divisions  needed  to  obtain  an  accurate  current  distribution  on  the 
dipole  could  be  reduced  by  assuming  a  more  complex  current  form  for  each 
dipole  section.  The  desirability  of  using  a  minimum  number  of  dipole  divisions 
is  two-fold.  The  most  obvious  reason  being  a  reduction  in  the  computer  time 
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B. _ --Continued 

required  for  a  given  degree  of  accuracy.  Although  a  savings  in  this  quantity 
is  not  particularly  important  for  the  single  dipole  solution,  it  is  essential  if  a 
multielement  structure  is  to  be  analyzed.  The  second  reason  for  requiring 
a  minimum  number  of  dipole  divisions,  although  less  obvious  than  the  first, 
is  equally  important.  In  reducing  the  theoretical  solution  presented  in  Section 
II  to  form  a  suitable  form  for  computer  use,  some  approximations  were  made  that 
require  the  length  of  each  dipole  division  to  be  much  greater  than  the  diameter 
of  the  dipole;  therefore,  there  is  a  maximum  number  of  dipole  divisions  that 
can  be  used  and  still  maintain  a  high  degree  of  accuracy  in  the  solution.  For 
thin  dipoles,  of  course,  the  problem  does  not  arise  since  an  accurate  solution 
can  be  obtained  by  using  a  number  of  dipole  divisions  far  below  the  maximum 
allowed. 
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Figures  III  - 1  and  IH-2  illustrate  the  current  distribution  on  the  half¬ 
wave  dipole  for  different  values  of  N  and  for  both  the  constant  current  model  on 
each  division  and  the  sinusoidal  current  model  on  each  division.  For  comparison, 
the  dipole  current  distribution  calculated  by  R.  W.P.  King*  is  also  shown  on  each 
graph.  The  advantage  of  using  the  sinusoidal  current  model  for  each  segment 
of  the  dipole  is  seen  by  comparing  the  two  figures.  With  the  sinusoidal  current 
form,  the  resulting  current  distributions  approach  that  obtained  by  King  with 
fewer  dipole  divisions  and  remain  nearer  King's  distribution  over  a  wider  range 
of  N  than  do  the  distributions  obtained  by  using  the  constant  current  form  for 

each  division.  It  is  also  interesting  to  note  that  the  maximum  current  is  not  at  the  center  of 
the  dipole  as  wculd  be  expected.  This  is  explained  by  the  fact  that  the  particular  dipole  used 
is  relatively  thick  and  thus  appears  electrically  longer  than  half  a  wavelength.  Current  dis¬ 
tributions  have  been  calculated  for  much  thinner  dipoles  and  the  current  maximum  does 
appear  at  the  center  of  the  dipole  as  expected.  (See  Figure  III  -  3 ) 

The  gain  of  this  antenna  was  calculated  using  the  method  described  in  ths  Sec  an  on 

2 

theory  and  was  2.15  db,  the  value  predicted  by  Jordan. 

C.  TWO  DIP6LES 

Having  verified  that  the  theory  was  correctly  applied  to  a  single  dipole,  the  next  step 
was  to  try  solving  the  current  distributions  on  two  dipoles  connected  by  a  transmission  line. 

Several  cases  of  two  dipoles  in  an  array  have  been  programmed  and  current  distributions 
found .  In  many  instances  much  useful  information  such  as  variations  in  input  impedance  with 
changing  dipole  spacings  and  lengths  was  found,  but  for  the  most  part  these  results  only 
served  as  a  check  on  the  theory  and  programming  procedure  and  were  not  investigated 
further. 

There  is  one  particular  case  of  general  interest,  however,  which  should  suffice  as  a 
representative  example  of  the  many  two-dipole  arrays  which  have  been  programmed.  The 
example  in  mind  is  an  array  of  two  half-wave  dipoles  separatee,  by  one-half  wavelength  and 
fed  in  phase  to  give  a  broadside  radiating  structure .  The  currents  on  both  dipoles  are 
identical.  The  dipoles  are  fed  by  a  single  voltage  generator  at  the  base  of  one  dipole,  with 
a  section  of  transmission  line  one-half  wavelength  long  connecting  the  remaining  dipole. 

To  obtain  "in  phase"  voltages  at  the  bas  e  of  each  dipole,  the  transmission  line  was 
switched  so  that  there  was  an  additional  180 -degree  phase  shift  at  the  second  dipole  besides 
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Figure  HI- 3.  Dipole  Current  Distribution  (Thin  dipole  -  fl  =  20) 
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the  180-degree  phase  shift  introduced  by  the  length  of  transmission  line.  The  gain  for 

3 

this  array  is  6.  029  db  and  is  consistent  with  values  found  by  using  other  methods.  The 
input  impedance  to  each  dipole  is  62  +  J6.  3  ohms  and  is  comparable  to  the  value  of 
60  +j4  obtained  by  Kraus. 


D.  LOG-PERIODIC  ANTENNA 

In  this  section  probably  the  most  important  results  of  the  study  are  presented. 

These  are  the  computed  values  of  current  distribution  and  radiated  power  for  an  entire 

log-periodic  antenna.  After  many  trial  computer  runs,  a  satisfactory  program  was 

finally  achieved  for  the  LP  antenna  and  the  results  of  this  program  were  compared  to 

data  measured  from  existing  antenna  models. 

The  model  that  was  analyzed  on  the  computer  was  a  10-element  switched-phase 

log-periodic  antenna  (bandwidth  250-500  MHz)  with  a  ~  45  and  T*  =  0.  853.  The  antenna 

feed  structure  was  a  137-ohm  two-wire  line.  Two  computer  runs  were  made  for  two 

different  frequencies.  Figures  III-4  and  III- 5  show  the  current  distribution  at  frequencies 

of  300  MHz  and  450  MHz,  respectively.  The  currents  are  nearly  as  would  be  expected 

by  antenna  engineers.  The  highest  currents  appear  on  the  elements  near  resonance  and 

on  the  next  smaller  elements;  i.  e.  ,  the  active  region  includes  the  resonant  element  plus 

the  next  several  smaller  elements.  The  phase  progression  from  one  element  to  ’he  next 

,  4 

is  also  in  agreement  with  measurements  made  previously. 


Figures  HI-6  through  IH-9  show  a  comparison  between  measured  and  calculated 
radiation  patterns  at  frequencies  of.  300  MHz  and  450  MHz  for  this  log -periodic 
dipole  array. 

Figures  IH-10  through  III-14  show  the  current  distribution  and  far-field 
patterns  (measured  and  calculated)  of  a  6-element  log-periodic  at  operating 
frequencies  of  650  and  850  MHz.  The  design  band  limits  are  550  and  1000  MHz. 

Figures  III- 15  through  IH-19  sho\v  the  current  distribution  and  far-field 
patterns  for  a  10- element  version  of  the  antenna  given  above.  In  comparing  the 
currents  shown  in  Figures  IH-10  and  HI-15  one  must  note  the  difference  in 
maximum  values.  The  major  reason  is  that  both  antennas  theoretically  were 
fed  with  ideal  1-volt  sources,  thus  aHowing  the  antenna  with  the  lower  input 
impedance  to  draw  more  current.  In  the  calculation  of  absolute  gain  as  described 
:»  Section  II,  this  effect  is  accounted  for  automatically. 
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Figure  III- 6.  Log-Periodic  Azimuth  Pattern  -  300  MHz 
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Figure  III-8.  Log-Periodic  Azimuth  Pattern  -  450  MHz 
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Figure  HI-10.  Current  Amplitude  and  Phase  Distribution 
6 -element  Log  Periodic.  T  =  0.  85  >  ct  =  90° 


FREQUENCY  =  850  MHz 
PATTERN  =  E  PLANE 
6  ELEMENTS 


180 


Figure  III- 12.  Log  Periodic  E-Plane  Pattern  -  850  MHz 
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Figure  III-13.Log  Periodic  H-Plane  Pattern  -  650  MHz 
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III- 14.  Log  Periodic  H-Plane  Pattern  -  850  MHz 
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Figure  III-17.  Log  Periodic  H-Plane  Pattern  -  b50  MHz 
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Figure  M-18.L,0g  Periodic  E-Plane  Pattern  -  850  MHz 


Figure  HI-19.  Log  Periodic  H-Plane  Pattern  -  850  MHz 


D.  --Continued. 

The  patterns  shown  in  Figures  III  -  6  -  9 ,  11-14,  and  16-19  are 
probably  the  most  conclusive  proof  available  to  show  that  the  theory  and  com¬ 
puter  program  are  valid  when  applied  to  a  log-periodic  antenna.  The  discrepan¬ 
cies  between  calculated  and  measured  data  that  show  up  in  the  back  lobes  of  the 
patterns  are  due  to  environmental  conditions  at  the  test  site  and  to  the  lack  of 
perfection  in  model  building  techniques.  One  argument  that  indicates  that  the 
physical  model  is  at  fault  instead  of  the  mathematical  model,  is  the  fact  that 
the  measured  data  is  not  symmetrical  where  the  differences  in  the  two  sets  of 
data  occur.  For  future  work,  higher  quality  models  will  be  used  and  tested  in 
a  controlled  environment. 

Figures  III-20  and  III-21  show  the  current  distribution  and  far- 
field  patterns  of  a  log  periodic  dipole  array  when  the  *  angle  is  45  degrees. 
Antennas  with  angles  this  high  are  sometimes  used  in  broadside  or  billboard 
arrays.  Since  the  computer  outputs  are  symmetrical  about  the  0-90  degree 
line,  only  half  of  the  E-  and  H-plane  patterns  are  included  in  Figure  HI-21  and 
many  of  the  succeeding  patterns. 

Figures  HI-22  through  III-24  show  the  current  distribution  and 
far  field  patterns  of  a  10-element  array  designed  for  maximum  gain  using  the 
techniques  developed  by  R.  L.  Carrel.  The  maximum  gain  was  found  to  be  10.  7 
dbi  while  Carrel's  method  predicted  a  gain  to  11  dbi. 

Figures  III- 25a  through  c  show  that  the  computer  program  that 

has  been  developed  will  predict  the  existence  of  the  so-called  3/2  moding.  Three- 
halves  moding  occurs  when  only  part  of  the  available  energy  is  radiated  in  the 
active  region  around  the  half  wave  element  with  the  remainder  being  radiated 
in  the  region  of  the  three-halves  wave  element.  The  effect  shown  probably 
would  have  been  more  pronounced  if  the  operating  frequency  had  been  a  little 
higher. 
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Figure  IH-23.  Log  Periodic  H -Plane  Patte 


Figure  III-25a  Current  Amplitude  and  Phase  Distribution  Log  Periodic  with  3^/2  Moding 


FREQUENCY  =  1950  MHz 
PATTERN  =  H-PLANE 


E. 


LOG -PERIODIC  ARRAYS 


As  a  final  step  in  the  computer  solution  of  log-periodic  antennas 
an  array  of  LP's  was  analyzed.  The  arrays  were  of  two  types  --  the  DuHamel 
array  where  all  log  periodic  antennas  making  up  the  array  are  identical  and  the 
Mei  array  where  all  elements  of  equal  length  lie  in  a  straight  line,  thereby 
making  it  necessary  for  each  LP  in  the  array  to  be  of  a  different  design  (See 
Figure  III-26).  The  DuHamel  array  has  a  circular  phase  front  whereas  the 
Mei  array  has  a  flat  phase  front.  Of  the  computer  runs  made,  the  DuHamel 
array  was  the  most  often  used  while  the  Mei  array  was  used  only  occasionally 
as  a  check  against  the  DuHamel  designs. 

Two  basic  LP  designs  were  used  as  elements  in  forming  the  LP 
arrays.  These  were  a  high  gain  LP  (10  dipoles,  X  =  6.  37°,  =  0.  96)  and  a 

medium  gain  LP  (9  dipoles,  =  44.  34°,  =0.  78).  Both  medium  and  high 

gain  antennas  were  designed  to  operate  in  the  upper  UHF  region  (length  of  longest 
dipole  =0.  25  meters)  so  models  could  be  constructed  and  tested  easily  if  a 
comparison  between  computer  and  empirical  data  was  desired.  This  choice  in 
design  frequencies  is  not  a  restriction  for  antennas  operating  in  free  space, 
however,  since  the  computer  program  in  this  mode  woi'ks  only  with  relative 
and  not  absolute  values,  i.  e.  ,  a  half  wavelength  antenna  looks  electrically  the 
same  at  all  frequencies.  The  only  case  where  actual  frequencies  must  be 
analyzed  in  the  computer  program  is  that  in  which  the  antenna  is  to  be  placed 
over  real  earth  which  is  frequency  sensitive  in  the  absolute  sense.  For  the 
sake  of  convenience,  all  patterns  for  the  medium  and  high  gain  antenna  arrays 
in  the  remainder  of  the  report  were  run  at  frequencies  of  1100  MHz  and 
665  MHz,  respectively.  See  Figures  HI-27  and  IH-28  for  single  LP  patterns  at 
these  frequencies. 
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In  switching  from  single  LP  antennas  to  LP  arrays  the  interests  in 
electrical  characteristics  changes  also.  Previously,  current  distributions 
were  of  prime  concern,  l  .t  now  the  gain  and  pattern  coverage  that  can  be 
achieved  with  different  antenna  designs  is  of  interest.  The  major  factors 
affecting  the  gain  and  pattern  coverage  of  an  array  are  the  individual  LP  designs 
of  which  we  have  chosen  two  as  being  representative,  the  manner  in  which  the 
LP's  are  arrayed  and  the  number  of  LP's  in  the  array.  The  particular  array, 
DuHamel  or  Mei,  once  chosen,  leaves  only  the  spacing  between  LP's  and  the 
number  of  LP's  in  the  array  as  the  variables  to  be  optimized. 

The  process  used  to  find  an  optimum  set  of  array  variables  was  one 
of  choosing  the  spacing  to  be  used  between  LP's  and  then  to  make  computer 
runs  with  varying  numbers  of  LP's  in  tne  array.  Figure  III-29  shows  a  graph 
of  maximum  gain  versus  the  number  of  LP's  in  the  array  for  various  DuHamel 
spacings.  The  term  LP  spacing  as  used  here  means  the  distance  in  wavelengths 
between  the  half  wave  dipoles  of  neighboring  LP's  in  the  array.  The  data  for 
this  graph  was  obtained  from  computer  runs  for  high  gain  LP's  in  an  array. 

It  is  interesting  to  note  that  the  array  gain  does  not  necessarily  increase  with 
an  increase  in  the  number  of  LP’s  in  the  array,  but  often  decreases.  If  one 
visualizes  the  geometry  of  an  LP  array  he  sees  that  each  additional-  element 
placed  ,n  the  array  must  be  positioned  at  an  ever  increasing  angle  from  the 
direction  of  radiation  for  the  array.  The  radiation  from  these  off-axis  LP's 
contributes  significantly  to  the  side  lobe  level  thereby  decreasing  the  amount  of 
gain  that  could  be  realized  with  the  additional  element.  This  increase  in  sidelobe 
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Figure  III-29  High  Gain  IP's  in  a  DuHamel  Array 


E. _ --Continued 

level  coupled  with  that  due  to  arraying  when  some  of  the  LP  elements  in  the 
array  are  separated  by  more  than  a  wavelength  can  cause  an  overall  decrease 
in  array  gain,  especially  when  the  LP  to  LP  spacing  in  the  array  becomes  large. 

Figure  M-30  is  a  graph,  similar  to  that  of  Figure  III- 29 ,  for  medium 
gain  LP's  in  an  array.  This  graph  shows  trends  identical  to  those  of  the  previous 
graph  with  the  only  difference  being  an  overall  lower  gain  for  the  medium  gain 
LP  arrays. 

The  two  graphs  just  mentioned  do  not  tell  the  complete  story  of  arraying 

trends;  however,  a  closer  look  at  some  radiation  patterns  indicates  that  a 

particular  array  may  have  more  gain  than  another,  but  its  radiation  coverage 

could  be  completely  unacceptable.  For  example.  Figure  III -29  shows  that  a 

DuHamel  array  with  four  LP's  spaced  1.  5  wavelengths  apart  has  a  higher 

maximum  gain  than  a  a  milar  array  with  elements  spaced  1.  25  wavelengths 

apart,  but  a  glance  at  the  radiation  patterns  shows  that  the  peak  of  the  beam 

(actually  two  peaks  --  one  on  either  side  of  boresight)  occurs  about  6  degrees 

off  boresight  for  the  higher  gain  array  (Figures  III-31a,  b).  See  Part  IV  - 

Appendix  for  a  complete  set  of  E-and  H-pIane  radiation  patterns  and  gains  of  all 

designs  run  on  the  computer.  Even  though  the  1.  5  wavelength  spaced  array  has 

higher  maximum  gain,  the  1.  25  wavelength  spaced  array  has  a  higher  boresight 
gain. 

In  all  previous  examples  the  arrays  analyzed  have  been  of  the 
DuHamel  or  circular  phase  front  type.  The  Mei  array,  however,  offers  an 
alternate  design  approach  which  seems  to  be  somewhat  superior  to  the  DuHamel 
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Figure  III-30  Medium  Gain  IP's  in  a  DuHamel  Array 
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array.  Figure  III-32  Is  a  bar  graph  showing  a  comparison  of  the  two  arrays 
for  different  designs  and  as  can  be  seen  the  Mei  array  is  significantly  higher 
in  gain  for  all  designs.  The  designs  used  for  this  comparison  were  not  preselected, 
but  were  chosen  simply  because  they  happened  to  be  the  only  designs  for  which 
both  array  types  were  run  on  the  computer.  The  comparison  is  by  no  means 
exhaustive  in  the  number  of  different  designs  used,  but  the  results  are  so  de¬ 
cidedly  one  way  that  the  Mei  array  could  safely  be  called  the  superior  array 
type  at  least  in  terms  of  array  gain.  One  also  notes  that  the  gain  of  the  Mei 
array  does  not  decrease  in  going  from  four  to  five  element  arrays  unlike  that  of 
the  DuHamel  designs.  There  is  nothing  to  indicate  that  this  trend  won't  reverse 
with  the  addition  of  the  sixth  element,  but  at  least  up  to  five  elements  the  gain 
of  the  Mei  array  increases  with  additional  elements.  See  Figure  III-33  for  the 
radiation  pattern  of  a  5  element  Mei  array-the  highest  gain  array  run. 

Before  concluding  this  section  on  the  comparison  of  array  types  we 
note  that  the  one  comparison  involving  medium  gain  LP's  shows  the  Mei  array 
with  higher  gain  than  the  DuHamel  array,  but  with  the  radiation  in  the  opposite 
direction.  Even  though  the  Mei  array  gain  is  higher,  the  fact  that  it  radiates  in 
the  backward  direction  seems  to  indicate  an  unstable  design.  Unfortunately 
there  are  no  other  computer  runs  on  medium  gain  LP  Mei  arrays  available  to 
indicate  whether  or  not  this  instability  is  a  property  of  the  array  type. 

In  our  final  step  in  the  analysis  of  LP  arrays  we  have  placed  the 
arrays  over  real  earth.  Due  to  an  oversight  on  the  part  of  the  author,  the 
computer  runs  were  made  at  UHF  instead  of  HF  and,  because  of  a  shortage 
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Figure  III-32  Comparison  of  DuHamel  and  Mei  Arrays 
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of  time,  reruns  on  these  arrays  at  the  proper  frequencies  were  not  obtained. 

The  two  cases  of  arrays  over  real  earth  that  were  run  are  of  some  interest, 

however,  and  are  shown  in  Figures  III-34a  and  b.  The  ground  chosen  for  the 

-14 

runs  had  a  conductivity,  W  ,  of  5  x  10  emu  and  a  relative  dielectric  constant, 
£  ,  of  15  -  a  rather  mediocre  earth.  The  gain  of  these  arrays  is  slightly  less 
than  there  counterparts  in  free  space  (See  Part  IV)  and  if  one  were  to  speculate 
he  could  draw  the  conclusion  that  this  decrease  in  gain  would  be  less  at  HF 
frequencies.  Speculating  further  we  can  say  that  most  probably  th*  tave-off 
angle  at  the  lower  frequencies  would  be  somewhat  higher  than  th'  .  i,  the  UHF 
case  at  hand,  although  there  is  no  way  of  specifying  exactly  what  it  would  be. 

The  parameters  affecting  the  radiation  characteristics  of  an  HF  array  are  so 
numerous  that  drawing  exact  conclusions  from  the  available  UHF  runs  would 
not  only  be  meaningless,  but  would  be  impossible.  It  is  only  possible  to 
speculate  about  the  direction  of  change  that  might  occur  in  going  from  UHF  to 
HF. 

Until  now  this  section  has  simply  provided  the  results  of  the  various 
computer  runs  on  LP  arrays  with  only  some  of  the  more  interesting  results  being 
pointed  out.  It  is  time  to  draw  some  conclusions  and  apply  these  conclusions 
to  a  specific  array  design.  At  the  beginning  of  this  study  a  set  of  design  goals 
to  be  achieved  was  specified.  With  the  limited  data  at  hand  we  can  see  which 
goals  can  be  fulfilled,  which  goals  seem  likely  to  be  fulfilled  and  which  goals 
will  take  considerably  more  study  to  determine  whether  or  not  they  are  possible 
to  meet. 
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One  of  the  foremost  design  goals  is  the  gain  which  has  been 
specified  as  a  minimum  of  17  dbi  at  the  beam  peak.  A  5  element  Mei 
array  with  3/4  wavelength  spacing  has  a  free  space  gain  of  16.  9  dbi,  nearly 
enough  to  meet  the  requirements.  However,  the  Mei  array  seems  to 
increase  in  gain  with  additional  elements  so  that  a  minimum  gain  of  17  dbi 
should  be  easily  achieved  with  a  6  or  7  element  array  even  over  real  earth. 

The  second  design  goal,  equally  important  as  gain,  is  the 
coverage  that  can  be  obtained.  The  azimuth  coverage  was  specified  as  a  16 
degree  beamwidth  at  and  near  the  peak  of  the  beam  in  elevation.  Since  the 
azimuth  coverage  should  change  only  slightly  between  a  free  space  situation 
and  over  real  earth,  the  free  space  patterns  available  should  give  a  good 
indication  as  to  the  coverage  that  can  be  achieved  by  an  HF  antenna.  Another 
look  at  the  same  5  element  Mei  array  shows  a  beamwidth  of  15  degrees  with  the 
first  sidelobe  17  db  below  the  beam  peak.  The  beamwidth  is  only  a  degree 
under  the  design  goal,  hardly  enough  to  bother  with  since  the  beamwidth  will 
probably  change  more  than  a  degree  with  a  change  in  frequency.  The  side- 
lobe  level  of  -17  db  is  2  db  better  than  the  -15  db  specified  in  the  design  goal 
requirements. 

The  elevation  coverage  is  to  have  a  beamwidth  from  2  to  16  degrees. 
We  can  see  from  the  UHF  pattern  that  this  goal  cannot  be  achieved  with  a  simple 
HF  antenna  above  ground,  there  must  ba  a  ground  screen  in  front  of  the  array 
to  lower  the  beam  peak.  No  more  can  be  said  about  the  elevation  pattern  or 
screen  effects  since  the  HF  information  is  not  available 
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By  using  the  data  that  is  available  on  the  five  element  Mei  array, 
we  can  now  proceed  with  the  design  of  an  optimum  3-30  MHz  HF  array.  As 
mentioned  previously,  the  number  of  elements  would  be  increased  from  five 
to  six  with  an  inter -element  spacing  of  0.  75  wavelengths. 

With  a  T  of  0.  96  each  LP  element  of  the  array  would  consist  of 
62  dipoles  for  a  total  of  372  dipoles  for  the  array.  The  longest  dipole  would 
be  164  feet  tall  while  the  shortest  would  be  13.  59  feet  tall.  Four  dipoles 
have  been  added  to  the  high  frequency  end  of  each  LP  element  to  insure 
that  maximum  gain  is  achieved  at  30  MHz. 

The  spacing  between  the  longest  dipoles  of  adjacent  LP  elements 
would  be  0.  75  X  or  246  feet  for  a  total  array  width  of  1230  feet  at  the  low 
frequency  end.  Similarly,  the  spacing  between  the  shortest  dipoles  of 
adjacent  LP  elements  would  be  20.4  feet  giving  an  array  width  of  102  feet 
at  the  high  frequency  end. 

Since  we  are  designing  a  Mei  array,  the  spacing  between  dipoles 
of  an  LP  element  is  dependent  on  element  placement  in  the  array  (See  Figure 
IH-26).  Although  the  dipole  spacings  differ  from  element  to  element,  only  the 
inter  dipole  spacings  of  an  imaginary  LP  in  the  center  of  the  array  need 
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to  be  specified  since  all  dipoles  of  the  same  length  on  the  various  LP 
elements  will  fall  on  a  line  perpendicular  to  this  imaginary  centerline. 

Of  course,  if  there  were  an  odd  number  of  LP  elements  in  the  array,  then 
there  would  be  an  actual  center  LP  element.  The  best  choice  for  the 
imaginary  center  LP  element  would  be  the  high  gain  LP  discussed  pre¬ 
viously  with  an  o^.  angle  of  6.  37°.  Using  this  information  and  some  simple 
geometry,  the  spacing  between  the  longest  and  the  next  longest  element  in 
the  imaginary  LP  can  be  found  from  the  expression 
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(153) 


where  S  is  the  desired  spacing  andjc*««vis  the  length  of  the  longest  element. 
For  the  particular  array  at  hand  this  spacing  is  59  feet.  Making  use  of 
geometrical  progression  formulas  one  can  determine  the  total  length,  L, 
of  the  array  from 
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where  is  the  length  of  the  smallest  element.  The  total  length  of 


our  array  is  1352.  1  feet. 


E. _ --Continued 

Until  now  only  the  dimensions  of  the  longest  dipole  and  the 
spacing  between  the  two  longest  dipoles  in  the  imaginary  center  LP  have 
been  presented.  Obtaining  the  dimensions  of  the  remainder  of  the  elements 
and  spacings  is  simply  a  matter  of  successive  multiplication  by  . 

The  only  remaining  parameter  that  must  be  specified  for  a 
complete  array  design,  is  the  line  length  between  adjacent  dipoles  in  an  LP 
element.  Although  the  physical  spacing  between  dipoles  of  an  LP  element 
depends  on  the  placement  of  the  element  in  the  array,  the  line  length  con¬ 
necting  any  two  dipoles  in  an  LP  element  must  be  the  same  for  similar  dipoles 
in  each  LP  element  in  the  array  to  maintain  a  flat  phase  front.  The  line 
lengths  can  be  obtained  by  letting  the  electrical  spacing  equal  the  physical 
spacing  for  the  outer  most  elements  in  the  array  since  the  outer  elements  will 
have  the  largest  separation  between  dipoles.  Again  using  some  simple  geometry 
we  find  that  the  line  length  from  the  largest  to  next  largest  dipole  in  the  outer 
LP  elements  is  63.  9  feet.  To  obtain  the  other  line  lengths  from  dipole  to 
dipole  in  the  LP's  all  that  is  required  is  successive  multiplication  of  the 
largest  line  length  by  ^  . 

The  design  of  the  HF  Mei  array  is  now  complete  and  may  be 
summarized  as  follows: 

Frequency  3-30  MHz 

No.  of  LP  Elements  6 

62 


No.  of  Dipoles/Element 
No.  of  Dipoles/Array 
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T  0.96 

Spacing  Between  Elements  246  ft.  Low  Frequency  End 

20.  4  ft.  High  Frequency  End 

Array  Width  1230  ft.  Low  Frequency  End 

102  ft.  High  Frequency  End 

Array  Length  1352.  1  ft. 

Largest  Dipole -to -Dipole  Spacing 

on  Centerline  59  ft. 

Largest  Line  Length  63.  9  ft. 

The  characteristic  impedance  of  the  transmission  line  was  60  ohms  for  the 

computer  runs  used  in  designing  this  array.  However,  no  attempt  has  been 

made  to  optimire  this  impedance,  therefore  no  transmission  line  impedance 

will  be  specified  in  this  design.  Similarly,  there  is  not  enough  information 

available  to  specify  the  type  and  size  ground  plane  necessary  to  make  this 

array  meet  the  elevation  coverage  specified. 

There  are  other  design  goals  such  as  VSWR  that  were  specified 
at  the  beginning  of  the  study,  but  the  program  has  changed  scope  to  such  an 
extent  that  these  goals  are  no  longer  in  the  realm  of  the  study.  Such  design 

goals  are  secondary,  however,  compared  to  the  primary  objective  of  proper 
gain  and  coverage. 
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Although  not  every  aspect  of  the  original  study  has  been  considered 
because  of  the  redirection  in  emphasis  from  antenna  physical  modeling  to 
mathematical  modeling,  it  is  felt  that  an  optimum  HF  array  design  has  been 
achieved  to  the  extent  intended  in  the  study.  During  the  course  of  the  study 
an  extremely  valuable  tool  in  the  form  of  a  computer  L.P  array  analysis  program 
has  been  developed  and  furnished  to  the  customer  for  optimizing  arrays  with 
different  design  goals. 

The  mathematical  model  developed  under  this  study  was  the  best 
at  the  time  of  its  writing,  but  use  of  the  program  and  an  analysis  of  the  results 
show  that  there  are  changes  that  would  improve  and  expand  this  model.  It  is 
felt  that  these  changes  would  greatly  benefit  the  user  of  such  a  model  and  that 
a  future  study  be  undertaken  to  incorporate  any  desireable  changes  in  the  model 
and  computer  program. 
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Figure  I V  A-l  High  Gain  Ii>  -  Gain  11.2  db 
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Figure  IV  b-5  2  Elenent  Array  -  1.5  Wavelength  Spacing  -  Gain  14.0  db 
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Figure  IV  B-15  5  Element  Array  -  1.5  Wavelength  Spacing  -  Gain  11.8  db 
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Figure  IV  G-3  2  Element  Array  -  0.5  Wavelength  Spaaing  -  Gain  9.4  db 
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Figure  IV  E-l  4  Element  Array  -  .5  Wavelength  Spacing  -  Gain  11. 3  db 
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the  antenna  all  electrical  characteristics  of  the  antenna,  such  as  radiation  patterns 
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space  and  real  earth  calculations  arc  presented  and  an  optimum  antenna  design  is 
chosen  from  this  information. 
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